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Preface 


At  the  request  of  the  Defense  Atomic  Support 
Agency,  Oak  Ridge  National  Laboratory  has 
undertaken  the  preparation  of  a  handbook  to  aid 
engineers  charged  with  the  responsibility  of 
designing  shields  to  protect  military  equipment 
and  personnel  in  the  vicinity  of  a  nuclear  weapons 
burst.  This  document  constitutes  one  chapter  of 
the  Handbook,  with  other  chapters  to  be  published 
as  they  are  completed.  These  will  include  a 
chapter  defining  the  radiation  sources  insofar  as 
is  possible  and  practicable  (Chapter  2),  a  chapter 
outlining  methods  for  calculating  the  attenuation 
of  weapons  radiation  through  various  media 
(Chapter  3),  and  a  chapter  presenting  radiation 
albedos  (Chapter  4).  These  four  chapters,  together 
with  an  introductory  first  chapter,  will  eventually 
be  bound  as  Volume  I  of  the  Handbook.  Volume  II 
will  consist  of  two  or  more  additional  chapters 
presenting  engineering  design  methods  that  are 
based  on  the  more  sophisticated  techniques 
described  in  Volume  I.  The  intent  is  that  the 
shield  designer  will  use  Volume  I  as  a  textbook 
and  ready  reference  and  Volume  II  as  a  guide  for 
handling  most  of  the  problems  with  which  he  will 
be  confronted. 

In  order  to  prepare  this  Handbook,  it  has  been 
necessary  for  Oak  Ridge  National  Laboratory  to 
obtain  the  assistance  of  several  consultants  and 
subcontractors.  For  this  chapter  on  ducts,  for 
example,  Wade  E.  Selph  of  Radiation  Research 
Associates,  Inc.,  performed  the  initial  literature 
search  and  prepared  the  draft  with  which  the 
co-author  and  the  editors  worked.  Other  chapters 
will  similarly  represent  a  cooperative  effort  of 
ORNL  and  other  organizations. 


As  is  always  the  case  for  handbooks,  the 
authors  and  editors  are  relying  heavily  on 
suggestions,  reviews,  and  criticisms  of  others  as 
an  aid  in  the  development  of  the  various  chapters. 
The  list  of  individuals  who  have  contributed  in 
this  manner  has  already  grown  very  large,  and  it 
would  be  almost  impossible  to  acknowledge  each 
person  here.  We  do,  however,  wish  to  express 
appreciation  to  Lt.  Col.  Charles  D.  Daniel  and 
Lt.  Col.  William  A.  Alfonte,  who  as  past  and 
present  DASA  Shielding  Project  Officers  have 
handled  the  administration  of  the  contract  and 
assisted  in  establishing  the  scope  of  the  Hand¬ 
book.  We  also  wish  to  acknowledge  the  assistance 
of  several  persons  at  Oak  Ridge  National  Labo¬ 
ratory  who  by  virtue  of  having  worked  in  the 
radiation  shielding  field  for  a  number  of  years  are 
able  to  serve  as  on-the-spot  authorities  to  help 
resolve  problem  areas  as  they-.arise  in  the  various 
chapters.  Those  who  have  acted  in  this  capacity 
for  the  duct  chapter  are  V.  R.  Cain,  R,  E.  Maerker, 
and  D.  K.  Trubey. 

Finally,  we  wish  to  thank  Mrs.  Virginia  M. 
Hamrick,  who  by  carefully  reading  each  draft  of 
this  chapter,  including  galley  and  page  proofs, 
has  both  improved  the  rhetoric  and  helped  eliminate 
some  of  the  usual  errors  that  always  creep  into 
formal  publications. 

December  1966 

The  Radiation  Shielding  Information  Center  grate¬ 
fully  acknowledges  the  interest  and  efforts  of  Lt. 
Col.  G.  C.  Reinhardt  and  Capt.  R.  W.  Enz  of  the 
Radiation  Physics  Branch  of  the  Defense  Atomic 
Support  Agency  for  making  it  possible  for  this 
work  to  be  reprinted  and  widely  distributed. 
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5.0 


Introduction 


Most  of  the  shields  that  are  designed  to  protect 
military  equipment  or  personnel  from  weapons 
radiations  will  contain  air-filled  holes  of  some 
type.  They  may  consist  of  access  ways  to  ac¬ 
commodate  ventilating  ducts  or  other  utility  pipes, 
passageways  to  allow  personnel  to  enter  the 
structure,  or  distributed  voids  resulting  from  the 
use  of  nonhomogeneous  material  in  the  structure. 
Since  radiation  traverses  air-filled  regions  es¬ 
sentially  unimpeded,  these  irregularities  in  an 
otherwise  adequate  shield  can  represent  a  serious 
problem  for  the  shield  designer. 

The  most  effective  method  for  reducing  the 
amount  of  radiation  that  travels  through  such 
openings  is  to  design  the  ducts  or  passageways 
so  that  they  do  not  penetrate  straight  through  the 
shield,  as,  for  example,  ducts  which  contain  one 
or  more  bends.  Radiation  traveling  through  ducts 
of  this  kind  collides  with  some  of  the  surrounding 
shield  material  and  consequently  is  attenuated  to 
some  small  fraction  of  the  incident  quantity.  Even 
if  during  such  collisions  particles  are  scattered 
back  into  the  duct,  their  angles  of  reflection  and 
their  energies  following  reflection  will  probably 
prevent  them  from  traveling  on  in  the  forward 
direction.  Thus,  the  introduction  of  a  bend  will 
eliminate  a  significant  fraction  of  the  radiation 
which  has  entered  the  mouth  of  the  duct.  This  is 
especially  true  for  gamma  rays,  since  high-energy 
gamma  rays  scattered  at  large  angles  suffer  a 
considerable  reduction  in  energy. 

When  the  radiation  being  considered  is  neutrons, 
the  problem  is  more  severe  in  that  neutrons  can 
undergo  many  scatterings  without  being  absorbed 
or  substantially  degraded  in  energy,  the  number  of 
scatterings  depending,  of  course,  on  the  compo¬ 
sition  of  the  shield  and  the  initial  energy  of  the 
neutron.  Consequently  neutrons  have  a  much  higher 
probability  for  “streaming”  through  shield  pene¬ 
trations  than  gamma  rays  do.  The  problem  is 
further  complicated  by  the  fact  that  neutrons  can 
introduce  new  sources  of  gamma  rays:  those 
emitted  when  neutrons  are  captured  and  those 
produced  when  neutrons  suffer  inelastic  scattering. 
These  neutron-induced  gamma  rays,  like  those 
from  the  original  source,  are,  of  course,  also  re¬ 
duced  in  intensity  by  the  introduction  of  bends. 


From  the  above  discussion  it  is  apparent  that  in 
a  practical  shield  pipes  should  follow  irregular 
paths  and  passageways  should  be  designed  with 
one  or  more  right-angle  turns.  Any  removable 
plugs  required  in  a  shield  should  be  of  the 
stepped-plug  type  so  that  radiation  streaming 
around  the  plug  will  strike  shielding  material 
before  it  penetrates  very  far  into  the  structure. 
But  even  with  all  these  precautions,  the  radiation 
transmitted  through  holes  may  easily  constitute 
the  major  fraction  of  the  dose  penetrating  a  shield. 
Thus  in  order  to  avoid  over-designing  and  thereby 
increasing  the  costs  associated  with  constructing 
ducts  and  passageways,  the  shield  designer  must 
be  able  to  predict  the  effect  of  various  types  of 
openings  for  a  given  set  of  conditions. 

Many  of  the  methods  that  have  been  developed 
for  estimating  the  effect  of  openings  in  reactor 
shields^ can  also  be  used  for  weapons  radiation 
shields.  (Initial  weapons  radiation  is  analogous 
to  radiation  from  an  operating  reactor,  and  fallout 
radiation  is  similar  to  that  from  a  shut-down  re¬ 
actor.)  Differences  will  exist  in  the  geometries 
and  energy  distributions  of  reactor  and  weapons 
sources  and,  in  general,  in  the  geometrical 
relationship  between  the  shield  and  the  openings 
in  the  shield,  but  they  can  be  accounted  for.  A 
weapons  radiation  shield,  as  is  demonstrated  by 
the  examples  given  in  this  chapter,  is  usually  con¬ 
sidered  to  be  a  thick-walled  structure  which  sur¬ 
rounds  a  central  cavity  and  is  penetrated  by  a 
duct  that  is  circular  or  rectangular  in  cross  sec¬ 
tion.  The  radiation  source  is  incident  on  the 
outside  of  the  structure  and  enters  the  structure 
both  by  penetrating  shielding  material  and  by 
following  the  path  of  the  duct. 

For  simple  ducts  (straight  ducts)  the  component 
that  will  probably  contribute  most  to  the  dose 
penetrating  the  structure  will  be  the  “line-of-sight” 
component  which  travels  uninhibited  from  the  mouth 
of  the  duct  to  the  central  cavity.  Analytic  methods 
for  calculating  this  component  are  presented  in 
Section  5.1.1.  Another  component  which  can  be 
important  for  simple  ducts  is  the  radiation  that 
reaches  the  duct  through  the  surrounding  shielding 
material,  either  as  uncollided  radiation  that  entered 
the  shield  in  the  direction  of  the  duct  or  as 
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radiation  that  is  scattered  by  the  shielding  ma¬ 
terial  into  the  duct.  These  components  are  gener¬ 
ally  handled  by  the  ray  analysis  technique  dis¬ 
cussed  in  Section  5.1.2. 

More  sophisticated  techniques  are  required  for 
treating  more  complicated  duct  geometries  or  for 
obtaining  higher  accuracy  for  the  simple  geometries. 
This  means  that  the  shield  designer  will  have  to 
resort  to  a  computer  program  to  perform  calculations 
that  involve  a  great  many  particle  scatterings  with 
consequent  changes  in  direction  and  energy 
degradation.  A  number  of  analog  Monte  Carlo 
computer  programs  are  available  for  such  calcula¬ 
tions,  but,  as  is  pointed  out  in  Section  5.1.3,  the 
amount  of  computer  time  required  to  use  them  is 
usually  prohibitive.  As  a  result,  modified  Monte 
Carlo  programs  have  been  developed  for  neutron 
transmission  calculations  which  do  not  attempt  to 
follow  particle  paths  through  the  shielding  material 


surrounding  the  duct  but  instead  employ  albedo 
theory  to  estimate  reflectance  from  the  duct  walls. 
These  programs  are  discussed  in  Section  5.1.4, 
along  with  other  albedo  techniques  that  have  been 
devised  for  estimating  gamma-ray  transmission. 

Finally,  this  chapter  presents  a  few  techniques 
for  calculating  the  effects  of  holes  or  depressions 
in  shields.  Both  single  large  voids  and  smaller 
voids  distributed  throughout  the  shield  are  con¬ 
sidered.  The  techniques  have  been  applied  almost 
exclusively  to  gamma  rays,  but  they  are  equally 
applicable  to  neutrons. 

The  applicability  of  most  of  the  methods 
presented  here  is  demonstrated  by  comparisons  of 
calculated  and  experimental  results.  In  a  few 
cases  additional  experimental  data  are  included 
to  provide  empirical  results  which  will  be  of  value 
for  future  reference. 
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5.1.  Methods  for  Calculating  Radiation 
Transmission  Through  Ducts 


5.1.1.  ANALYTIC  METHODS  FOR 
LINE-OF-SIGHT  RADIATION 

For  the  case  of  a  duct  that  penetrates  straight 
through  a  thick  shield  and  has  dimensions  that  are 
large  compared  with  the  mean  free  path  of  the 
penetrating  radiations,  the  “line-of-sight”  radia¬ 
tion,  that  is,  the  radiation  which  travels  directly 
from  the  source  to  the  detector  through  the  duct, 
can  be  the  most  significant  portion  of  the  total 
radiation  reaching  the  detector.  A  simple  example 
of  this  situation  is  that  of  a  disk  source  and  an 
adjacent  shield  which  has  a  cylindrical  duct 
passing  completely  through  it.  For  a  detector 
located  on  the  duct  axis  three  or  more  source 
diameters  from  the  source,  the  ratio  of  the  line-of- 
sight  radiation  through  the  duct  to  the  uncollided 
radiation  reaching  the  detector  through  the  shield 
is  approximately  the  ratio  of  the  cross-sectional 


area  of  the  duct  divided  by  the  area  of  the  disk 
source  times  the  attenuation  of  the  shield.  Thus 
if  the  area  of  the  source  is  ten  times  the  area  of 
the  duct  and  the  shield  is  approximately  three 
mean  free  paths  thick,  the  line-of-sight  component 
is  two-thirds  of  the  total  uncollided  radiation  reach¬ 
ing  the  detector. 

In  more  complex  situations,  for  example,  in  the 
case  of  a  weapon  burst  and  the  subsequent  scat¬ 
tering  of  the  radiation  before  it  enters  the  duct, 
the  source  can  be  treated  as  a  plane  source  at  a 
given  location  and  will  have  a  particular  angular 
distribution  dependent  upon  the  geometry  of  the 
actual  source  and  its  surrounding  medium.  If  it  is 
assumed  that  the  source  is  a  plane  source  at  the 
duct  opening  and  that  the  ratio  of  the  duct  length 
to  diameter  is  large,  all  line-of-sight  points  are 
approximately  equidistant  from  the  other  end  of 
the  duct  and  may  be  treated  collectively  as  an 
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equivalent  point  source.  If  the  plane  source  is 
treated  as  having  isotropic  current,  the  equivalent 
point  intensity  is  the  product  SN^,  where  is  the 
particle  emission  rate  per  unit  area  and  time  and  S 
is  the  area  on  the  source  plane  that  is  viewed 
through  the  duct  from  the  point  of  interest.  For 
nonisotropic  emitters  each  differential  area  ele¬ 
ment  must  be  weighted  by  the  angular  distribution 
in  performing  the  integral  over  the  source  area. 
For  complex  source  distributions  and  geometries, 
the  integral  over  the  source  must  be  evaluated 
numerically.  However,  for  many  geometries  and 
source  angular  distributions  which  may  be  ex¬ 
pressed  analytically,  it  is  possible  to  derive 
formulas  for  the  line-of-sight  flux  or  current  through 
an  opening. 

The  sections  below  present  several  methods  that 
have  been  developed  for  computing  the  line-of-sight 
flux  in  simple  ducts  of  various  geometries.  For 
these  cases  the  source  is  assumed  to  be  a  plane 
which  is  normal  to  the  axis  of  the  duct  and  covers 
the  duct  entrance.  Two  types  of  source  angular 
distributions  are  considered;  an  isotropic  distribu¬ 
tion  and  cosine  distributions.  Here  isotropic  dis¬ 
tribution  refers  to  radiation  that  is  assumed  to  be 
given  off  in  all  directions  in  the  forward  hemisphere 
with  equal  probability  but  with  no  emission  in  the 
backward  hemisphere.  Thus  if  a  surface  emission 
is  particles /cm ^  sec,  then  the  number  of  par¬ 
ticles  emitted  into  any  unit  solid  angle  is  N 
particles/cm  ^  steradian  sec. 

For  cosine  distributions  the  radiation  leakage 
through  a  unit  area  on  the  surface  of  a  self-absorb¬ 
ing  volume-distributed  source  is  most  intense  in 
the  direction  normal  to  the  surface,  and  the  varia¬ 
tion  with  the  angle  0  from  the  normal  is  approxi¬ 
mated  by  a  cosine  function,  in  which  case  the 
number  of  particles  emitted  per  unit  solid  angle  is 
cos  d/TT  particles/cm  ^  steradian  sec. 

Other  functions  have  been  used  to  describe 
sources  which  have  more  forwardly  peaked  angular 
distributions  than  that  given  by  cos  6.  In  general, 
most  of  these  can  be  represented  by  cos"  6,  where 
n  may  be  as  large  as  20.  The  cos"  6  function  has 
been  useful  in  representing  a  leakage  current  or  a 
flux  at  a  shield  surface  from  an  absorbing  source 
region  in  which  the  activity  per  unit  volume  in¬ 
creases  with  depth. 

Two  types  of  plane  sources  have  been  used  in 
analytical  expressions  derived  for  the  transmis¬ 
sion  of  uncollided  particles  through  ducts.  In  one 


type  the  source  strength  is  defined  by  a  particle 
emission  rate,  which  in  this  chapter  will  always 
be  designated  as  N ^  (particles/cm  ^  sec).  In  the 
other  type,  the  source  strength  is  defined  by  a 
flux,  which  will  be  designated  here  as  (par¬ 
ticles/cm  ^  sec).  This  second  source  is  the  more 
natural  one  for  use  in  shield  design.  It  is  used 
almost  exclusively  for  shelter  design  since  an 
analysis  of  the  transmission  of  radiation  from  a 
weapon  burst  through  ambient  atmosphere  will 
yield  a  flux  which  will  then  be  used  as  the  source 
at  the  duct  entrance.  (An  example  of  an  isotropic 
plane  source  described  by  a  particle  emission  rate 
is  a  converter  plate  such  as  the  one  that  was  used 
in  the  Lid  Tank  Shielding  Facility  at  Oak  Ridge 
National  Laboratory  as  a  source  for  shielding 
experiments.) 

In  calculating  the  flux  in  a  duct,  care  must  be 
taken  to  ensure  that  the  proper  source  angular  dis¬ 
tribution  is  used  in  a  given  situation  since  there 
is  a  basic  difference  in  the  two  sources.  This  is 
because  a  particle  emission  rate  is  a  current  and 
therefore  is  referenced  to  a  square  centimeter  of 
the  source  area,  whereas  a  flux  is  always  refer¬ 
enced  to  a  square  centimeter  of  area  that  is  normal 
to  the  particle  direction.  Therefore,  to  correct  for 
the  different  orientation  of  the  flux-source  refer¬ 
ence  area,  it  is  necessary  to  multiply  the  flux 
passing  through  that  area  by  cos  6  in  order  to 
project  the  reference  area  back  to  the  source  plane. 
That  is,  the  relationship  between  the  angular  cur¬ 
rent  /(cos  6)  from  the  first  type  of  source  and  the 
angular  flux  <I)(cos  &)  from  the  second  type  of 
source  is 


/(cos  0)  =  cos  6  <I>(cos  6)  .  (5.1) 

The  effect  of  the  difference  in  the  two  sources 
can  be  shown  by  the  following  example.  The  un¬ 
collided  flux  in  a  duct  is  to  be  calculated  for  a 
given  angular  distribution  g(cos  6)  (particles/cm^ 
steradian  sec)  specified  on  the  source  plane.  If 
g(cos  6)  is  the  normalized  angular  current  [the 
normalization  condition  on  g(cos  6)  is  g(cos  6) 
2tt  d(cos  6)  =  1],  then  the  flux  at  any  point  in  the 
duct  is 


I 


Ngg(cos  0) 


dS 


(5.2) 
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where  is  the  total  emission  rate,  dS  is  a  dif¬ 
ferential  area  on  the  source  plane  S,  and  r  is  the 
distance  between  dS  and  the  detector.  If,  on  the 
other  hand,  g(cos  0)  is  the  normalized  angular 
flux,  then  the  flux  in  the  duct  is 

cos  9 

<I>Q  g(cos  6)  ■■  dS  ,  (5.3) 

where  is  the  total  or  scalar  flux  at  the  source 
plane. 

A  general  relation  between  the  emission  rate  and 
the  flux  may  be  derived  for  the  very  useful  case  of 
the  cos"  6  distribution  by  defining  the  angular 
emission  rate  (angular  current)  as 


crossing  a  plane  S  (cross-sectional  area  of  duct 
normal  to  axis)  is  assumed  to  be  known  and  is 
designated  by  an  arbitrary  function  g(cos  6),  which 
is  the  angular  flux  and  represents  the  number  of 
particles  traveling  in  the  direction  9  with  respect 
to  the  duct  axis  per  unit  solid  angle,  unit  time,  and 
unit  area  normal  to  direction  9  (neutrons/cm  ^  sec 
steradian).  The  flux  at  any  point  along  the  axis  is 
expressed  by 

(D  =  g(cos  9)  d^Ji9)  ,  (5.7) 

where  d 0,^(9)  is  the  differential  solid  angle  sub¬ 
tended  from  an  isotropic  detector  by  the  surface 
dS  (see  Fig.  5.1)  and  is  equal  to 


/(cos  9)  = 


n  +  1 
277 


Ag  cos"  9  . 


(5.4) 


cos  9  dS 


That  Eq.  5.4  is  properly  normalized  is  seen  by 
calculating  the  total  emission  rate: 

cos"  9  277  d{cos  0)  =  .  (5.5) 

277  0  V  /  0  V  7 

The  corresponding  expression  for  the  total  flux, 
at  the  source  plane  is  then 

^  n  +  1  cos" 9 

- - N. -  277  d(cos  9) 

277  0  cos  (9 


n  +  1 


It  is  evident  that  Eq.  5.6  does  not  hold  when  n  =  0; 
in  fact,  the  flux  is  not  defined  on  the  source  plane 
when  the  emission  rate  is  isotropic. 

For  all  the  different  duct  geometries  included  in 
the  following  discussion,  the  flux  0  or  leakage 
current  /  at  the  exit  plane  of  the  duct  is  expressed 
in  the  same  units  as  the  source  strength,  e.g., 
particles/cm  ^  sec.  In  order  to  obtain  the  total 
leakage  flux  or  current  through  the  exit  plane  of 
the  duct,  it  is  necessary  to  integrate  <E>  or  /  over 
the  exit  area. 

Rectangular  Ducts.  —  The  work  of  Hubbell  et  ah* 
on  radiation  from  uniform  rectangular  sources  with 
arbitrary  angular  distributions  is  directly  applicable 
to  the  calculation  of  the  uncollided  flux  along  the 
axis  of  a  rectangular  duct.  In  that  work,  radiation 


A  solution  to  Eq.  5.7  was  developed  by  Hubbell 
et  al.  in  terms  of  completely  separable  source  and 
geometry  functions  and  was  expressed  as 

.  (5.8) 

/=0  ^ 
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Fig.  5.1.  Schematic  of  Source  Plane  and  Detector 
Position  for  a  Rectangular  Duct. 
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where  a  is  the  ratio  of  the  half  height  (H/2)  to  the 
length  (Z)  of  the  duct,  b  is  the  ratio  of  the  half 
width  (W'/2)  to  the  length,  and  gj  and  pj  are  coef¬ 
ficients  of  Legendre  polynomial  expansions  of  the 
source  and  geometry  functions  respectively.  With 
these  definitions  of  a  and  b,  Eq.  5.8  gives  the  flux 
at  a  corner  for  the  quarter  plane  limned  in  Fig.  5.1. 
The  flux  at  the  center  line  for  the  full  plane  is  ob¬ 
tained  by  multiplying  Eq.  5.8  by  4. 

The  coefficients  are  calculated  by 


and 


g(cos  6)  Pj(cos  6)  d(cos  0)  (5.9) 


2)  when  uniform  over  the  duct  mouth  are  given  be¬ 
low.  The  flux  at  a  corner  of  a  rectangle  with 
dimensions  W/2  x  H/2  is  one-fourth  of  that  cal¬ 
culated  by  Eqs.  5.13  through  5.15. 

For  n  =  0: 


<D 


2$ 


-“tan-' 


ab 


n 


^/n~a^~+~b^ 


(5.13) 


For  n  =  1: 

2$„ 


TT 


a  b 

— . -  tan  -  ' - 


+  a^ 


p/a,b)  =  /  P/cos  6)  dQJi6)  ,  (5.10) 

where  the  Pj(cos  0)’s  are  the  Legendre  polynomials. 

Although  the  solution  of  Eq.  5.10  is  in  terms  of  a 
finite  series,  the  terms  become  progressively  more 
complicated  with  increasing  /  and  only  the  first 
four  pj’s  were  evaluated  analytically.  Using  an 
electronic  computer,  Hubbell  et  al.  numerically 
evaluated  p^’s  for  0  ^  =  13  over  a  grid  of  a  and  b 

values  of  0.1  ^  6  =  a  ^  20.  Using  these  tabulated 
values  and  the  g^’s  determined  for  any  arbitrary 
^(cos  6),  the  flux  along  the  axis  or  at  the  comer 
of  a  duct  can  be  calculated  with  Eq.  5.8. 

For  the  special  and  useful  case  of  the  cos"  6 
distribution. 


b  a 

,  ■„  tan-' - - 

VTTP 


(5.14) 


For  n  =  2: 


$  = - 2. 


77 


ab 


+  \  1  -I-  1  +  6^ 


.  (5.15) 


As  the  distance  between  the  source  and  detector 
becomes  large  with  respect  to  cross-section  dimen¬ 
sions,  the  inverse  tangent  in  Eq.  5.13  approaches 
ab  and  the  terms  in  the  brackets  in  Eqs.  5.14  and 
5.15  approach  2ab  or  WH/2Z^.  In  general,  it  can 
be  shown  that  for  long  thin  ducts  of  rectangular 
cross  section 


g(cos  0)  = 


(n  +  1)  ^0 

277 


cos"  0  , 


(5.11) 


(n  +  1) 
277Z2 


(5.16) 


where  <t)^  is  the  source  strength  that  is  represented 
by  a  scalar  flux  obtained  by  integrating  the  angular 
flux  over  all  angles  in  the  direction  of  the  duct 
entrance.  Consequently,  Eq.  5.11  is  properly 
normalized  since 

'  n  +  1 

-  cos"  0  277  d(cos  0)  =  1  .  (5.12) 

277 

For  integral  values  of  n,  Eq.  5.11  can  be  ex¬ 
panded  into  a  finite  number  of  Legendre  poly¬ 
nomials,  and,  in  principle,  an  analytical  solution 
to  Eq.  5.8  can  be  obtained.  However,  the  solution 
rapidly  becomes  unwieldly  for  progressively  hi^er 
values  of  n.  The  fluxes  along  the  axis  of  a  rec¬ 
tangular  duct  for  the  isotropic  flux  source  (n  =  0) 
and  for  more  forwardly  peaked  sources  (n  =  1  and 


Equation  5.16  represents  a  lower  limit  for  a  point 
source  (emitting  only  into  the  forward  hemisphere) 
and  also  for  a  plane  source  (for  a  detector  far 
away  from  the  duct  entrance)  since  the  flux  scat¬ 
tered  from  the  wall  of  the  duct  is  not  included. 
In  many  practical  cases  the  contribution  to  the 
dose  by  the  scattered  flux  will  be  of  the  same 
order  as  or  smaller  than  the  contribution  by  the 
uncollided  flux,  (The  dose  is  calculated  by  sum¬ 
ming  the  product  of  the  flux  and  the  appropriate 
dose  response  function  for  each  energy  group  of 
neutrons.) 

Rectangular  Slots.  —  The  geometry  used  in  the 
derivations  of  the  approximate  equations  for  a 
rectangular  slot  is  shown  in  Fig.  5.2.  The  dimen¬ 
sions  of  the  duct  adjacent  to  the  source  plane  are 
W  hy  H  and  the  distance  from  the  source  plane  to 
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Fig.  5.2.  Rectangular  Slot  Geometry.  Source  plane 
and  exit  plane  are  normal  to  Z  direction. 

the  exit  end  of  the  slot  is  Z.  When  H  »  W  (say, 
H/W  >  5),  the  exposed  source  can  be  approximated 
by  a  line  source,  and  the  line-of-sight  flux  0  and 
the  leakage  current  J  at  any  point  P  along  the 
horizontal  axis  on  the  exit  end  of  the  slot  by 


If  the  slot  dimensions  are  such  that  Z  »  H  and 

0j  +  02  the  above  equations  become 

N.W 

<D  .  -  °  , 

'  2Z 

(5.21) 

“  ttZ  ’ 

(5.22) 

2N^W 

$  -  “  , 

=  ttZ 

(5.23) 

J  ^0^ 

2Z 

(5.24) 

The  rectangular  slot  can  be  considered  a  special 
case  of  the  rectangular  duct  for  the  condition  when 
W  -r  0,  and  Eq.  5.19  can  be  derived  by  taking  the 
limit  of  Eq.  5.13  as  W  ^  0.  In  the  limiting  process 
an  apparent  anomaly  occurs;  the  answer  is  in  terms 
of  $  ,  the  flux  at  the  source.  However,  the  fluxes 
on  a  line  source  for  both  isotropic  and  cosine 
distributions  are  always  infinite.  Consequently, 
in  taking  the  limit  of  Eq.  5.13  as  W  -♦  0,  should 
be  converted  to  an  emission  rate  i2N^). 

Cylindrical  Ducts.  -  Consider  a  cylindrical  duct 
of  radius  a  normal  and  adjacent  to  a  plane  source 
with  the  detector  a  distance  Z  from  the  source 
plane  and  on  the  duct  axis  (Fig.  5.3).  For  an  iso¬ 
tropic  current  or  emission  rate  on  the  source  plane. 
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N.W 

N  W 

] .  =  — ^  (sin  0,  +  sin  0  )  ; 
'  27tZ  ^  ^ 


(5.17) 

(5.18) 


$  = — ^  (sin  6.  +  sin  0  )  ,  (5.19) 


N  W 

J  = — 5—  [sin  20,  +  sin  20  +  2(0  +  02)] ;  (5.20) 

47rZ  '  ^  ^ 


where  N ^  is  the  source  emission  rate  per  unit  area 
or  current,  0  is  expressed  in  radians,  and  the  sub¬ 
scripts  i  and  c  refer  to  isotropic  and  cosine  sources 
respectively. 
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the  line-of-sight  flux  at  any  point  P  is  given  by 
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3>.  =  ^  In  [1  +  (a/Z)2]  ,  (5.25) 

and  the  current  parallel  to  the  axis  of  the  cylinder 
is  given  by 


1  - 


Vl  +  (a/Z)2 

For  a  cosine  current  source 

=  2^0  '  ^ 


(5.26) 


Vl  +(a/zy 


(5.27) 


and 


Z2[l  +  (a/Z)2]  • 


(5.28) 


Fig.  5.4.  Geometry  for  Cylindrical  Annulus  Adjacent 
to  Plane  Source. 


axis  may  be  approximated  by 


Equations  5.27  and  5.28  can  be  expressed  in  terms 
of  an  isotropic  flux  at  the  source  plane  by  letting 
2^0  = 

For  Z  »  a,  say,  Z/a  >  5,  it  can  be  shown  that 
these  equations  approach 


(5.31) 


$  =  7  = 

i  •'  I 


Ng(7Ta^) 

2ttZ^ 


2Z2 


(5.29) 


which  is  the  same  result  for  a  point  source  of 
strength  N ^(jra^),  and 


N  (TTa^) 

//V-  O'  0 

r* 


■7TZ2 


(5.30) 


If  the  condition  of  a  long  duct  of  small  cross- 
sectional  area  is  met,  Eqs.  5.29  and  5.30  will 
also  hold  for  a  straight  duct  of  any  cross-sectional 
geometry  if  rra^  is  replaced  by  the  cross-sectional 
area  of  the  duct. 

Long  Cylindrical  Annulus.  —  The  geometry  for  a 
long  cylindrical  annular  duct  is  shown  in  Fig.  5.4. 
The  duct  is  bounded  by  two  cylindrical  surfaces 
having  radii  and  a^.  For  the  conditions  Z  » 
—  a^  and  32  »  a,  -  a^  and  an  isotropic  source 
N the  average  line-of-sight  flux  $  along  the  duct 


The  average  leakage  current  J  for  the  isotropic 
source  can  also  be  obtained  from  Eq.  5.31  by 
setting  M  =  1.  The  equation^  for  the  flux  at  any 
radius  r  (a^  =  r  ^  a^)  within  the  annulus  is 

MN  a^ 

<I)(r,Z)  = - ^  {1  -  (a  /a  )2  cos"!  (a  /r) 

ttZ  ^ 

+  cos-*  (Sj/a^) 

-(a/a,)[l_(a/a,)2]i/2S,  (5.32) 

where  M  =  1  for  an  isotropic  source  and  M  =  2  for  a 
cosine  source,  and  r  and  Z  are  the  cylindrical 
coordinates  of  the  point  P  at  which  the  flux  is 
given.  Equation  5.32  is  equivalent  to 

MN^S 

0(r,Z)  = - 2^  ,  (5.33) 

27rZ2 

where  S  is  the  area  on  the  source  plane  that  is 
viewed  from  the  point  P(r,Z). 
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5.1.2.  RAY  ANALYSIS  TECHNIQUE 

General  Description.  —  In  addition  to  the  line- 
of-sight  radiation,  the  flux  reaching  a  detector 
located  near  the  exit  end  of  a  simple  duct  includes 
radiation  that  enters  the  duct  through  its  walls.* 
This  contribution  consists  of  two  components: 
radiation  that  travels  directly  from  the  source  to 
the  detector  through  shielding  material  without  an 
interaction  (uncollided  flux)  and  radiation  that 
scatters  in  the  direction  of  the  detector  as  a  re¬ 
sult  of  interactions  with  the  shield.  A  method 
useful  in  calculating  these  components  is  the  ray 
analysis  method,  which  has  frequently  been 
applied  to  shield  attenuation  calculations. 

The  basic  assumption  of  the  ray  ana'ysis  tech¬ 
nique,  also  referred  to  as  the  point  kernel  method, 
is  that  the  radiation  transmission  is  a  function 
only  of  the  path  lengths  through  each  material  or 
void  encountered  along  a  straight  line  between  the 
source  point  and  the  detector.  This  assumption  is 
entirely  correct  for  predicting  the  uncollided  flux, 
and  for  geometries  in  which  the  scattered  flux  be¬ 
comes  important,  the  collided  fluxes  can  be  ac¬ 
counted  for  to  a  first  approximation  by  using 
buildup  factors  for  gamma  rays  and  removal  cross 
sections  for  neutrons.  Although  the  accuracy  of 
the  method  is  highly  dependent  on  the  particular 
configuration  being  considered,  it  can  be  applied 
to  simple  ducts  with  considerable  confidence 
since  the  radiation  reaching  the  detector  largely 
consists  of  the  uncollided  component. 

To  illustrate  the  method,  consider,  for  example, 
the  case  of  a  point  detector  P  removed  a  distance 
r  from  a  point  source  at  Q,  where  the  straight-line 
path  between  the  two  points  passes  through  both 
shielding  material  and  void  (see  Fig.  5.5).  The 
flux  at  the  detector  is  given  by 

,  (5.34) 

^  47TT 

where  N  is  the  point-source  strength  and  K  is  the 
material  attenuation  kernel  for  all  the  materials 
located  between  the  source  and  the  dose  point. 


*Radiation  that  enters  the  mouth  of  the  duct  in  a 
direction  other  than  the  direction  of  the  detector  but 
subsequently  scatters  in  the  duct  wall  toward  the  de¬ 
tector  also  contributes  to  the  flux.  For  simple  ducts 
the  contribution  from  this  source  is  usually  small,  but 
for  ducts  with  bends  it  becomes  the  dominant  compo¬ 
nent  as  is  discussed  in  Sections  5.1.3  and  5.1.4. 
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Fig.  5.5.  Geometry  for  Ray  Analysis  Calculation  of 
Radiation  Component  Arriving  Through  Wall  of  Duct. 

If  only  the  uncollided  flux  at  P  is  considered. 


where  I.  is  the  mean  free  path  of  a  particle  and  t. 
is  the  straight-line  path  through  the  ith  material. 
When  the  radiation  being  considered  consists  of 
neutrons,  l/I-  is  usually  expressed  in  terms  of  the 
total  macroscopic  cross  section  2^,  (cm“^)  and  t. 

is  given  in  centimeters.  When  gamma  rays  are  being 
treated,  l/L  is  expressed  as  the  mass  absorption 
coefficient  (cm^/g)  and  t.  is  in  units  of 

cm^/g.  If  the  thickness  t.  is  expressed  in 
centimeters,  then  1/1  is  replaced  by  the  linear 
attenuation  coefficient  (cm~^),  which  is  the  total 
macroscopic  cross  section  for  gamma  rays. 

Equation  5.35  can  be  modified  to  include  a  first- 
order  approximation  of  the  scattered  gamma 
radiation  by  multiplying  the  right-hand  side  by  a 
buildup  factor: 

K  .  .  .)  exp  L  t.  ^  .  (5.36) 

The  equation  can  be  modified  similarly  to  include 
a  first-order  approximation  of  scattered  neutrons 
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by  substituting  the  neutron  removal  cross  section* 
S^forS,: 

K  =  exp  -  I  2^  t.  .  (5.37) 

In  the  above  discussion  the  assumption  is  made 
that  the  flux  at  the  point  P  is  the  same  as  it  would 
be  if  all  the  particles  traveling  from  Q  to  P,  in¬ 
cluding  those  that  are  scattered,  encounter  ef¬ 
fectively  the  same  amount  of  attenuating  material 
that  exists  along  the  line  from  Q  to  P.  If 
significantly  less  material  is  encountered  by  some 
particles,  the  total  flux  at  P  determined  by  Eq. 
5.34  may  be  underestimated.  Conversely,  if 
significantly  more  material  is  encountered,  the 
flux  may  be  overestimated.  Also,  Eq.  5.34  neglects 
the  effects  of  material-void  interfaces  or  interfaces 
between  dissimilar  materials  encountered  along  the 
ray  and  in  the  surrounding  regions. 

In  the  use  of  the  ray  analysis  method  to  integrate 
over  a  source  plane  or  source  volume,  the  possible 
inaccuracy  involved  in  the  contribution  of  an 
individual  source  point  may  be  masked  by  the 
contributions  from  other  points.  There  are,  of 
course,  complex  geometries  and  material  combi¬ 
nations  for  which  the  simple  ray  analysis  tech¬ 
niques  do  not  predict  the  radiation  fluxes  with 
any  satisfactory  degree  of  accuracy.  However,  ray 
analysis  is  in  many  cases  quite  adequate  and  is 
an  easily  applied  method. 

There  are  also  cases  in  which  the  accuracy  re¬ 
quirements  and  the  complexity  of  the  problem 
geometry  may  prohibit  the  use  of  the  approximations 
given  in  this  section  except  as  first  estimates. 
In  such  situations  numerical  integration  may  be 
used  to  obtain  a  ray  analysis  solution  either  by 
hand  or  by  machine  calculations.  The  decision 
to  use  a  computer  will,  in  general,  be  based  on  the 
number  of  calculations  that  must  be  performed  to 
give  the  desired  degree  of  accuracy.  Several 
computer  programs  have  been  developed  which 
facilitate  point-kernel  integrations.  Some  of  their 
common  features,  as  well  as  brief  descriptions  of 
a  few  specific  programs  that  are  available,  are 
given  in  the  appendix  to  this  chapter. 

Application  to  Cylindrical  Ducts.  —  For  the  case 
in  which  the  duct  shown  in  Fig. 5.5  is  a  cylindrical 


^Removal  cross  sections  are  valid  only  for  those 
cases  in  which  the  shielding  materials  are  followed  by 
a  thick  layer  of  water  or,  to  a  lesser  extent,  for  those 
cases  in  which  mixtures  containing  hydrogenous  ma¬ 
terials  are  used  (see  Chapter  3). 


duct  of  radius  a  penetrating  an  infinite  slab  shield 
adjacent  to  an  infinite  plane  isotropic  source,  the 
uncollided  flux  at  P  that  arrives  through  the  duct 
walls  is 

oo  ,  dr 

f  ,  (5.38) 

r  0  Jr  r 

0 

where 

=  surface  emission  rate  or  current, 

1(t)  =  path  length  in  the  shield  which  lies  along  r 


fi  =  linear  attenuation  coefficient. 

The  corresponding  flux  for  a  cosine  source  is 
given  by 

oo  ,  Z 

$c=2Ao/  — dr  .  (5.39) 

Equations  5.38  and  5.39  are  not  amenable  to 
straightforward  analytical  solutions,  although 
values  of  the  integrals  may  be  found  in  tables  of 
mathematical  functions.  Certain  approximations, 
based  on  a  Taylor  expansion  of  /(r),  have  been 
given  by  Chase,*  but  they  are  valid  only  for  large 
values  of  a/Z.  Trubey®  evaluated  the  equations 
numerically  for  a  unit  surface  source  intensity 
(for  =  1).  The  results  are  shown  in  Tables  5.1 
and  5.2  as  a  function  of  the  ratio  of  the  duct  radius 
to  the  length  (a/Z)  and  the  shield  thickness  /rZ, 
where  pZ  is  measured  in  relaxation  lengths  (or 
mean  free  paths  determined  from  material  cross 
sections). 

The  data  for  the  uncollided  flux  given  in  Tables 
5.1  and  5.2  may  be  used  to  estimate  the  total 
radiation  arriving  through  the  duct  wall,  provided 
that  the  parameters  a/Z  and  pZ  are  measured  in 
terms  of  an  effective  relaxation  length,  A  =  1/p, 
which  accounts  for  the  radiation  scattered  in  the 
direction  of  P  by  the  shield.  Relaxation  length 
values  which  include  scattering  may  be  obtained 
from  experimental  data  or  from  basic  calculations 
of  penetration  through  the  material  by  Monte  Carlo 
or  moments  methods.  When  such  data  are  used 
they  should  be  for  the  same  material  and  the  same 
energy  (or  spectrum)  as  those  being  considered  in 
the  duct  penetration  problem.  Neutron  removal 
cross  sections  or  gamma-ray  buildup  factors,  both 
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Table  5.1.  Uncollided  Flux  at  Exit  of  Cylindrical  Duct  Due 
to  Radiation  Arriving  Through  Duct  Walls  (Isotropic  Source)^ 


Flux  Per  Unit  Surface  Source  Intensity  for  Shield  Thickness  of 

\z) 

flZ  =  0.1* 

flZ  =  0.2 

liZ=  0.5 

HZ=  1.0 

/lZ=  2.0 

fiZ  =  5.0 

0.001 

1.823 

1.223 

0.5602 

0.2198 

0.04912 

1.166  X  10~^ 

0.002 

1.823 

1.223 

0.5607 

0.2207 

0.04935 

1.184  X  10“^ 

0.005 

1.824 

1.224 

0.5620 

0.2215 

0.05002 

1.242  X  10“^ 

0.01 

1.825 

1.226 

0.5643 

0.2235 

0.05119 

1.349  X  10“^ 

0.02 

1.828 

1.229 

0.5688 

0.2278 

0.05363 

1.614  X  10“® 

0.05 

1.834 

1.239 

0.5820 

0.2406 

0.06171 

2.977  X  lO"^ 

0.1 

1.842 

1.254 

0.6024 

0.2613 

0.07621 

7.042  X  10“^ 

0.2 

1.851 

1.272 

0.6360 

0.2982 

0.1049 

1.795  X  10“^ 

0.5 

1.825 

1.277 

0.6841 

0.3659 

0.1687 

5.223  X  10“^ 

0.75 

1.765 

1.240 

0.6804 

0.3819 

0.1923 

6.990  X  10“^ 

1.0 

1.689 

1.185 

0.6560 

0.3767 

0.1980 

7.753  X  10~^ 

®From  D.  K.  Trubey,  A  Calculation  of  Radiation  Penetration  of  Cylindrical  Duct  Walls,  ORNL-CF-63-2- 
64  (Feb.  28,  1963). 

^’Number  of  relaxation  lengths. 

Table  5.2.  Uncollided  Flux  ot  Exit  of  Cylindrical 
to  Radiation  Arriving  Through  Duct  Walls  (Cosine 

Duct  Due 

Source)® 

Flux  Per  Unit  Surface  Source  Intensity  for  Shield  Thickness  of 

(z) 

^lZ  =  0.1* 

llZ=:  0.2 

fiZ=  0.5 

ijtZ  =  1.0 

IJ.Z=2.0 

llZ=5.0 

0.001 

1.445 

1.149 

0.6539 

0.2976 

0.07544 

0.002026 

0.002 

1.446 

1.149 

0.6545 

0.2983 

0.07582 

0.002059 

0.005 

1.446 

1.150 

0.6564 

0.3002 

0.07698 

0.002166 

0.01 

1.447 

1.152 

0.6595 

0.3034 

0.07896 

0.002365 

0.02 

1.450 

1.156 

0.6657 

0.3100 

0.08314 

0.002862 

0.05 

1.455 

1.166 

0.6836 

0.3298 

0.09703 

0.005458 

0.1 

1.459 

1.179 

0.7095 

0.3610 

0.1217 

0.01323 

0.2 

1.453 

1.188 

0.7462 

0.4120 

0.1682 

0.03343 

0.5 

1.339 

1.116 

0.7510 

0.4714 

0.2494 

0.08728 

0.75 

1.190 

0.9949 

0.6827 

0.4465 

0.2553 

0.1036 

1.0 

1.036 

0.8643 

0.5953 

0.3956 

0.2340 

0.1014 

^From  D.  K.  Trubey,  A  Calculation  of  Radiation  Penetration  of  Cylindrical  Duct  Walls,  ORNL-CF-63-2' 
64  (Feb.  28,  1963). 

^Number  of  relaxation  lengths. 
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Fig»  5.6»  Effect  of  Duct  Radius  and  Shield  Thickness 
on  Total  Uncollided  Flux  Reaching  Exit  of  Cylindrical 
Duct  (Isotropic  Source). 

of  which  are  derived  from  basic  penetration  data, 
may  be  used  in  obtaining  values  of  A.  In  this  use 
of  the  data,  the  assumption  is  made  that  the  total 
attenuation  is  exponential;  that  is,  the  relaxation 
length  does  not  vary  with  penetration  depth. 

Figures  5.6  and  5.7,  which  are  based  on  Trubey’s 
data,®  show  the  variation  with  shield  thickness 
and  duct  size  of  the  total  flux,  including  the  line- 
of-sight  component,  reaching  the  center  of  the  duct 
exit.  All  parameters  are  measured  in  terms  of  the 
shield  material  relaxation  length.  The  flux  is 
given  for  an  infinite  plane  source  of  unit  source 
strength  (e.g.,  1  particle/cm^  sec). 

The  relative  importance  of  the  line-of-sight 
radiation  discussed  in  Section  5.1.1  and  of  the 
uncollided  radiation  that  leaks  through  a  duct  wall 
is  illustrated  by  some  additional  calculations  by 
Trubey  given  in  Table  5.3.  Here  it  is  apparent 
that  the  ratio  of  the  line-of-sight  component  to  the 
total  uncollided  flux  depends  strongly  on  the  ratio 
of  the  duct  radius  to  the  length  (a/Z)  and  on  the 
shield  thickness  /zZ.  These  ratios  represent  an 
upper  limit  since  the  buildup  of  the  flux  due  to 
scattering  has  not  been  considered.  The  ratios 
are  applicable  to  either  fast  neutrons  or  gamma 
rays  provided  the  shield  thickness  is  defined  in 
terms  of  the  number  of  relaxation  lengths.  An 
examination  of  the  data  in  Table  5.3  reveals  that 
the  radiation  which  penetrates  the  wall  of  a  duct 
dominates  over  that  from  the  source  in  the  mouth  of 


Fig.  5.7.  Effect  of  Duct  Radius  and  Shield  Thickness 
on  Total  Uncollided  Flux  Reaching  Exit  of  Cylindrical 
Duct  (Cosine  Source). 

the  duct  for  small-diameter  ducts  that  penetrate 
thin  shields.  As  the  duct  diameter  increases 
and/or  the  shield  thickness  increases,  the  line- 
of-sight  component  becomes  of  greater  importance. 

If  the  duct  is  filled  with  a  lightly  attenuating 
medium,  formulas  similar  to  those  in  preceding 
sections  could  be  derived  which  include  expo¬ 
nential  attenuation  through  the  duct  material  in  the 
kernel  integration.  This  conditiori  would  tend  to 
decrease  further  the  importance  of  the  source  at 
the  mouth  of  the  duct. 

Application  to  Partially  Penetrating  Cylindrical 
Ducts.  —  The  radiation  transmitted  through  the 
walls  of  a  cylindrical  duct  that  penetrates  only 
part  of  the  shield  and  is  not  adjacent  to  the  source 
may  be  calculated  by  the  ray  analysis  technique  in 
much  the  same  manner  as  for  a  completely 
penetrating  duct.  Consider,  for  example,  the  duct 
shown  in  Fig.  5.8a.  Equations  5.38  and  5.39  may 
be  applied  to  this  geometry  merely  by  changing  the 
lower  limit  of  integration  from  to  Z.  As  before, 
/(r)  is  that  portion  of  r  which  lies  inside  the  shield. 

An  alternative  to  numerical  integration  of  these 
equations  for  each  case  would  be  to  determine  the 
flux  penetrating  an  infinite  slab  of  thickness  t  (see 
Section  5.2.1),  use  this  as  a  source  strength  input 
to  the  ducted  shield  of  thickness  Z  —  t,  and  eval¬ 
uate  using  the  data  presented  in  the  previous  sec¬ 
tion  (Figs.  5.6  and  5.7).  In  the  absence  of  more 
definitive  data,  a  more  forwardly  peaked  distribution 
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Table  5.3. 

Ratio  of  Line-of-Sight  Radiotion  to 

Totol  Uncolli 

ided  Flux  in  Duct^ 

Fraction  o 

f  Total  Flux  for  Shield  Thickness  of 

(z) 

pZ=  0.1* 

jlZ=Q.2 

pZ  =  0.5 

yiZ  =  1.0 

pZ  =  2.0 

j^lZ=  5.0 

Isotropic 

Infinite-Plane 

Source 

0.001 

<10“® 

<10~® 

<10“® 

<10“® 

<10“'' 

<10“^ 

0.002 

<10“® 

<10“® 

<10“® 

<10“® 

<10“'' 

0.002 

0.005 

<10~® 

<10“'' 

<10“'' 

<10“'' 

<10“® 

0.010 

0.01 

<10“'^ 

<10“'' 

<10“'' 

10“^ 

0.001 

0.038 

0.02 

<10“^ 

<10“® 

10“3 

10“^ 

0.004 

0.11 

0.05 

<10“^ 

0.001 

0.002 

0.005 

0.020 

0.29 

0.1 

0.003 

0.004 

0.008 

0.019 

0.062 

0.42 

0.2 

0.011 

0.015 

0.031 

0.062 

0.16 

0.52 

0.5 

0.059 

0.083 

0.14 

0.23 

0.40 

0.68 

0.75 

0.11 

0.15 

0.25 

0.37 

0.53 

0.76 

1.0 

0.17 

0.23 

0.34 

0.48 

0.64 

0.82 

Cosine  1 

Infinite-Plane 

Source 

0.001 

<10~® 

<10“® 

<10“® 

<10“® 

A 

o 

1 

<10“^ 

0.002 

<10“® 

<10“® 

<10“® 

<10“'' 

A 

O 

1 

0.002 

0.005 

<10““* 

<10“'' 

<10“'' 

<10“'' 

A 

o 

1 

o.> 

0.011 

0.01 

<10“'' 

A 

o 

1 

<10“^ 

<10“^ 

0.001 

0.040 

0.02 

<10“^ 

<10“^ 

<10“^ 

0.002 

0.005 

0.12 

0.05 

0.002 

0.002 

0.004 

0.008 

0.025 

0.27 

0.1 

0.007 

0.008 

0.014 

0.027 

0.076 

0.43 

0.2 

0.026 

0.031 

0.050 

0.12 

0.19 

0.54 

0.5 

0.14 

0.16 

0.22 

0.31 

0.45 

0.71 

0.75 

0.25 

0.28 

0.37 

0.48 

0.61 

0.79 

1.0 

0.36 

0.40 

0.50 

0.60 

0.71 

0.85 

^From  D.  K.  Trubey,  A  Calculation  of  Radiation  Penetration  ot  Cylindrical  Duct  Walls,  ORNL-CF-63-2- 
64  (Feb.  28,  1963). 

^’Number  of  relaxation  lengths. 


such  as  a  cosine  distribution  should  be  used  at  the 
artificial  interface  if  t  is  greater  than  one  relaxa¬ 
tion  length.  Even  when  the  original  source  is 
isotropic,  the  radiation  will  have  taken  on  a  more 
directional  character  after  penetrating  a  relaxation 
length. 

Simpler  equations  may  be  derived  for  the  flux 
through  partially  penetrating  ducts  if  the  assump¬ 
tion  is  made  that  only  the  radiation  penetrating 


the  end  of  the  partial  duct  need  be  considered.  In 
this  case  the  flux  atP  in  Fig.  5.8a  may  be  estimated 
by 

$  -  E^iiit  sec  'I')]  ,  (5.40) 

where  *1*  =  tan~^[a/(Z  -  ()]•  Similarly,  the  flux  at 
P  due  to  radiation  entering  the  mouth  of  a  duct  that 
penetrates  only  part  of  the  shield  but  is  adjacent 
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Fig*  5.8*  Geometries  for  Roy  Analysis  Colcuiotions  of  Radiation  Component  Arriving  Through  Wails  of  Ducts  That 
Partially  Penetrate  Shield* 


to  the  source  (see  Fig*  5*8i)  is  given  by 

4)  sec 'P')]  .  (5.41) 

where  W  ^  =tan"'(a/Z)* 

For  the  same  shield  and  duct  dimensions,  the 
angle  ^ '  is  smaller  than  *P  and  the  peak  flux  in¬ 
crease  at  P  due  to  the  presence  of  the  duct  is 
greater  in  case  b  than  in  case  a*  For  the  approxi¬ 
mation  fit  «2Z^,  the  integrated  emergent  cur- 
rentmaybeapp  roxim  a  tedby7Ta^Af^e~'“^[a  ^  /2(Z  —  f)^] 
for  either  case  a  or  b.  Thus,  while  the  peak  current 
tends  to  be  greater  in  case  b,  the  integrated  cur¬ 
rent  increase  due  to  the  void  is  approximately  the 
same  in  both  cases* 

Comparison  with  Experiment*  —  The  ray  analysis 
technique  has  been  used  successfully  many  times 
to  predict  experimental  results  for  both  gamma-ray 
and  neutron  transmission  in  solid  unpenetrated 
shields*  In  complex  geometries  involving  ducts, 
the  technique  should  give  better  results  for  gamma 
rays  than  for  neutrons,  since  neutrons  on  the  average 
undergo  more  scatterings  before  they  are  absorbed 
or  escape*  However,  as  is  demonstrated  by  the 


comparisons  discussed  below,  the  technique  yields 
good  results  for  the  case  of  neutrons  transmitted 
through  straight  ducts  abutting  a  reactor* 

The  ray  analysis  technique  was  employed  by 
Benenson  and  Fasano®  to  analyze  experiments  at 
the  Brookhaven  National  Laboratory  Shielding 
Facility  on  the  transmission  of  fast  neutrons  through 
straight  cylindrical  ducts  in  a  water  shield*  Neu¬ 
trons  from  a  fission  source  plate  in  the  shielding 
facility  entered  the  bottom  of  a  water  tank  in  which 
a  duct  had  been  vertically  positioned  and  neutrons 
which  reached  the  detector  at  the  opposite  end  of 
the  duct  arrived  there  either  through  the  base  of  the 
duct  or  through  its  walls  after  passing  through  the 
surrounding  water*  The  ducts  were  2,  4,  6,  8,  12, 
and  16  in*  in  diameter*  The  ®^S(n,p)®^P  reaction 
was  used  to  measure  the  fast-neutron  flux,  and  an 
experimentally  determined  relaxation  length  was 
used  to  analyze  the  penetration  of  neutrons  above 
the  threshold  from  the  fission  source*  Figure 
5.9,  which  is  typical  of  the  results  of  this  work, 
shows  that  the  ray  analysis  method  apparently 
yields  good  results  when  applied  to  streaming  of 
very  penetrating  radiation  from  a  diffuse  source 
through  simple  ducts* 
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Fig.  5.9.  Comparison  o^  Ray  Analysis  Calculations 
with  Measurements  of  Fast-Neutron  Flux  Along  Axis  of 
8-ln.-diam  Straight  Cylindrical  Duct  In  a  Water  Shield: 
Fission  Source  at  Duct  Mouth.  (From  Benenson  and 
Fasano,  ref.  6.) 

Piercey  and  Bendall^  also  used  the  ray  analysis 
technique  to  calculate  the  flux  of  fast  neutrons 
transmitted  through  straight  cylindrical  ducts. 
Their  calculations  corresponded  to  measurements 
made  for  1-,  2-,  and  3^^-in.-diam  aluminum  ducts 
positioned  in  the  water  of  the  LIDO  Shielding 
Facility  at  Harwell,  England.  The  duct  lengths 
were  up  to  200  duct  radii.  As  in  the  experiment 
discussed  above,  the  fast  flux  was  determined  by 
the  ^^S(n,p)^^P  reaction.  The  calculations  pre¬ 
dicted  the  absolute  sulfur  reaction  rate  to  within 
a  factor  of  2  over  the  experimental  values.  The 
poorest  agreement  occurred  in  the  region  of  the 
duct  nearest  the  reactor.  It  was  concluded  that 
most  of  the  discrepancy  resulted  because  of  the 
inhomogeneity  of  the  system  and  particularly  be¬ 
cause  of  the  aluminum  layer  between  the  duct  face 
and  the  reactor.  Figure  5.10  shows  a  comparison 
of  the  calculated  and  experimental  results  as  an 
example  of  this  work. 


10‘*  10°  10'  10^ 

DISTANCE  FROM  DUCT  MOUTH  (IN  DUCT  RADII) 


Fig.  5.10.  Comparison  of  Ray  Analysis  Calculations 
with  Measurements  of  Fast-Neutron  Flux  Along  Axis  of 
3^-in.-diam  Straight  Cylindrical  Duct  in  a  Water  Shield: 
Duct  Mouth  Adjacent  to  Reactor  Core.  (From  Piercey 
and  Bendall,  ref.  7.) 


5.1.3.  ANALOG  MONTE  CARLO  METHOD 

General  Description.  —  It  is  shown  in  the  pre¬ 
ceding  sections  that  for  straight  ducts  in  simple 
geometry  the  most  important  components  contribut¬ 
ing  to  the  flux  at  the  detector  are  the  line-of-sight 
radiation  and  the  radiation  which  enters  the  duct 
through  its  wall  at  the  proper  angle  to  reach  the 
detector.  For  more  complicated  geometries  involv¬ 
ing  ducts  with  bends,  the  line-of-sight  component 
disappears  and  the  dominant  component  becomes 
the  radiation  that  is  transmitted  through  the  duct 
by  successive  scatterings  from  the  surrounding 
shield  material.  Calculations  of  this  contribution 
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require  more  sophisticated  techniques  than  those 
described  in  the  preceding  discussion,  since  on 
the  average  many  scatterings  are  involved  and  sev¬ 
eral  boundaries  are  crossed. 

A  method  that  is  potentially  exact  for  calculating 
the  radiation  transmitted  through  a  duct  is  the 
Monte  Carlo  method,  often  called  “an  experiment 
run  on  a  computer.”  In  general,  the  Monte  Carlo 
method  traces  the  path  of  a  particle  until,  based  on 
a  probability  distribution  along  the  path,  a  colli¬ 
sion  occurs.  The  nature  of  this  event  is  deter¬ 
mined  by  selecting  from  a  probability  distribution 
of  all  possible  interactions.  If  the  event  is  an  ab¬ 
sorption,  the  history  is  generally  terminated.  If  it 
is  a  scatter,  direction  and  distance  to  the  next 
collision  are  determined  and  the  process  is  re¬ 
peated  until  the  particle  is  absorbed  or  until  it  ex¬ 
ceeds  the  geometric  or  energy  bounds  of  the  pro¬ 
blem.  In  most  Monte  Carlo  programs  there  are  wide 
variations  in  the  statistical  techniques,  allowable 
events,  and  source-shield  geometry  which  may  be 
considered. 

There  ate  two  types  of  Monte  Carlo  methods 
that  can  be  applied  to  duct  transmission  problems; 
analog  Monte  Carlo  methods,  which  ate  the  methods 
referred  to  in  this  section,  and  albedo  Monte  Carlo 
methods,  which  are  described  in  Section  5.4.4. 
The  term  “analog”  is  used  to  designate  calcula¬ 
tions  in  which  the  model  solved  is  an  analog  of 
the  interactions  occurring  as  the  particles  traverse 
the  duct  and  confining  walls,  as  opposed  to  the 
albedo  calculations  in  which  the  particles  travers¬ 
ing  a  duct  are  followed  by  a  random-walk  process 
and  wall  interactions  are  represented  by  a  reflec¬ 
tion  coefficient. 

The  analog  Monte  Carlo  method  is  not  always 
practical  with  the  present  programs  and  comput¬ 
ing  machinery  because  of  the  machine  computing 
time  required;  nevertheless,  several  machine  pro¬ 
grams  that  have  been  developed  are  applicable 
to  this  type  of  calculation.  Included  among  these 
are  the  ADONIS,  L05,  COHORT,  and  05R  codes, 
all  of  which  are  described  briefly  in  the  appendix 
to  this  chapter. 

Comparison  with  Experiment.  —  Collins  and 
McCleary®  used  the  L05  Monte  Carlo  code  to  cal¬ 
culate  the  transmission  of  neutrons  from  a  Po-Be 
source  (3.08  x  10^  neutrons/sec)  through  12-in. - 
long,  3-  and  6-in.-diam  straight  cylindrical  ducts 
in  water,  with  the  detector  positioned  3  in.  from 
the  duct  mouth.  The  agreement  that  they  obtained 


Fig.  5.11.  Comparison  of  LOS  Analog  Monte  Carlo 
Calculations  with  Measurements  of  Fast-Neutron  Dose 
Rates  Along  Axes  of  12-in. -long  3-  and  6-in.-diam  Straight 
Cylindrical  Ducts  in  Water:  Po-Be  Neutron  Source. 

(From  Collins  and  McCleary,  ref.  8.) 

between  the  measured  and  calculated  data,  shown 
in  Fig.  5.11,  lends  confidence  in  the  ability  of 
the  L05  procedure  to  predict  neutron  fluxes  and 
dose  rates  beyond  a  straight  cylindrical  duct  in 
a  shield. 

The  L05  code  was  also  used  by  Marshall®  for 
calculations  corresponding  to  measurements  made 
of  radiation  transmitted  through  cylindrical  ducts 
penetrating  water  shields.  The  ducts  were  con¬ 
structed  of  3-in.-diam  aluminum  tubing  and  had 
one  bend  midway  along  their  lengths,  with  the 
bend  angle  varying  between  30  and  60°.  The 
sources  were  14-MeV  neutrons  and  ®°Co  gamma 
rays  positioned  on  the  duct  center  line  3  in.  from 
the  mouth.  The  detector  at  the  opposite  end  of 
the  duct  was  also  3  in.  from  the  duct  mouth. 

Figures  5.12  and  5.13  show  a  comparison  of 
some  typical  experimental  data  and  calculated 
results  for  the  case  of  a  30  °  bent  duct  through 
an  18-in.  water  shield.  In  general,  Marshall  found 
that  for  shielding  thicknesses  which  exceeded 
18  in.  the  calculated  dose  rates  were  in  good 
agreement  with  the  experimental  data  in  shape 
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Fig.  5.12.  Comparison  of  LOS  Analog  Monte  Carlo 
Calculations  with  Measurements  of  Fast-Neutron  Dose 
Rate  Transmitted  Through  30°  Bent  Duct  (3-in.  diameter 
by  18  in.)  in  Water  Shield;  14-MeV  Neutron  Source. 
(From  Marshall,  ref.  9.) 

but  were  too  low  in  magnitude.  He  attributed  this 
to  three  possible  causes:  the  number  of  histories 
used  in  the  calculations  was  insufficient,  there 
was  structure  surrounding  the  experimental  con¬ 
figuration  that  was  not  considered  in  the  calcula¬ 
tions,  and  there  were  errors  in  the  cross-section 
values  used. 

Gardner  and  Mettler^®  used  the  ADONIS  code 
to  predict  the  transmission  of  neutrons  from  a 
T((f,n)  reaction  source  (14.7  MeV)  through  11-in.- 
and  3-ft-square  concrete  ducts.  They  compared 
their  results  with  measurements  made  by  Doty 
for  an  essentially  isotropic  source  and  with  the 
albedo-type  calculations  made  by  Song^^  and  by 
Maerker.^'*  (Further  discussion  of  the  measure¬ 
ments  and  albedo  calculations  is  give  in  Section 
5.1.4.)  Maerker’s  calculations  were  for  lower 
source  energies  and  there  is  some  confusion  in 
regard  to  the  normalization  of  the  various  results, 
but  since  the  correlations  and  results  are  pre¬ 
liminary  in  nature,  Gardner  and  Mettler  concluded 


Fig.  5.13.  Comporrson  of  LOS  Analog  Monte  Carlo 
Calculations  with  Measurements  of  Gamma-Ray  Dose 
Rate  Transmitted  Through  30°  Bent  Duct  (3-in.  diameter 
by  18  in.)  in  Water  Shield:  ^®Co  Gamma-Ray  Source. 
(From  Marshall,  ref.  9.) 

that  reasonably  good  agreement  exists  between 
the  calculations  and  measurements.  They  also 
pointed  out  that  the  experimental  measurements 
were  subject  to  large  error  as  a  result  of  sta¬ 
tistical  variation  of  the  measured  dose  rates,  cal¬ 
ibration  of  the  instrument,  and  determination  of 
the  incident  dose,  and  claimed  only  an  order  of 
magnitude  accuracy. 

Gardner  and  Mettler  also  used  the  ADONIS  code 
to  study  the  effect  of  the  “corner  lip”  in  multi¬ 
legged  rectangular  ducts  (see  Section  5.1.4). 

5.1.4.  ALBEDO  METHODS 

General  Description.  —  As  is  pointed  out  in 
Section  5.1.3,  the  primary  means  by  which  radiation 
is  propagated  through  a  duct  with  bends  is  by  suc¬ 
cessive  scatterings  in  the  duct  walls.  The  most 
practical  way  to  evaluate  this  component  is  by 
using  an  albedo  method.  In  general,  this  method 
assumes  that  the  point  at  which  radiation  is  inci¬ 
dent  on  the  duct  wall  is  a  point  on  an  infinite 
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interface  and  that  a  certain  fraction  of  the  inci¬ 
dent  radiation  will  emerge  from  the  wall  at  the 
same  point.  This  fraction  is  called  the  albedo. 
In  duct  transmission  calculations  both  total  and 
differential  albedos  have  been  used.  For  best 
accuracy  the  albedos  must  be  differential  both  in 
energy  and  in  the  angles  of  incidence  and  reflec¬ 
tion.* 

The  simplest  example  of  an  albedo  type  of  duct 
transmission  calculation  is  one  which  considers 
radiation  that  has  been  reflected  from  a  duct  wall 
only  once  before  reaching  the  detector.  The  cal¬ 
culation  is  performed  by  integrating  a  “reflection 
kernel”  over  the  duct  surfaces  which  have  an  unob¬ 
structed  line  of  sight  to  both  the  source  and  the 
detector.  Consider,  for  example,  the  case  of  a 
point  source  and  simple  duct  shown  in  Fig.  5.14. 
The  dose  at  P  due  to  first-order  reflection  from  the 
wall  is 
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Fig.  5.14.  Geometry  for  Albedo  Calculation  of  Singly 
Reflected  Radiation  in  a  Simple  Duct  (Point  Source). 
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wall 

surface 


NpR(Eg)  cos  dS 
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(5.42) 


where 

Np  =  source  strength. 

Eg  =  source  energy, 

P(Eq)  =  flux-to-dose  conversion  factor, 

~ dose  albedo  (re¬ 
flected  current  per  incident  cur¬ 
rent),** 

cf)  =  azimuthal  angle  separating  the 
incident  and  reflected  rays, 

and  fj,  r^,  6^,  and  6  are  defined  in  Fig.  5.14.  The 
reflected  dose  must,  of  course,  be  added  to  the 
line-of-sight  dose  to  obtain  the  total  dose. 

The  Simon-Clifford  equations'®  for  neutrons  re¬ 
flected  through  ducts  having  cylindrical  geometry 
(discussed  later  in  this  section)  resulted  from  an 
integration  of  Eq.  5.42  for  an  assumed  analytical 
fit  to  CL.  For  cases  in  which  a  fit  to  their  ex¬ 
pression  for  a  is  possible,  the  Simon-Clifford  equa¬ 
tions  should  be  used.  Other  fits  may  be  derived 


which  would  allow  the  integration  to  be  performed 
for  simple  geometries. 

If  either  the  nature  of  a  or  the  geometry  prohibits 
direct  integration,  Eq.  5.42  may  be  solved  numeri¬ 
cally  by  dividing  the  duct  wall  into  small  scatter¬ 
ing  areas  and  assuming  that  all  points  within  an 
elemental  scattering  area  are  the  same  distance 
from  the  source  or  detector  and  that  they  have  the 
same  angles  of  incidence  and  reflection.  A  numeri¬ 
cal  evaluation  performed  by  LeDoux  and  Chilton'® 
led  to  a  set  of  equations  for  calculating  singly 
reflected  gamma  rays  propagated  through  bent  rec¬ 
tangular  ducts.  (The  LeDoux-Chilton  technique  is 
discussed  below.) 

In  the  case  of  a  broad-beam  source  incident  on 
the  duct  at  an  angle,  such  as  that  shown  in  Fig. 
5.15,  no  iine-of-sight  radiation  reaches  the  detec¬ 
tor  and  it  is  more  important  that  all  reflected  radi¬ 
ation  be  evaluated,  including  that  which  has  been 
reflected  more  than  once.  The  singly  reflected 
component  at  P  is  given  by 

Dp  =  j  ^0 

over  ^ 

illuminated 
area 


*Definitions  of  albedos  and  corresponding  albedo 
data  are  presented  in  Chapter  4  of  this  Handbook. 

**See  the  definitions  given  in  Chapter  4;  various  other 
albedo  forms  may  also  be  used. 


(5.43) 

where  D.  is  the  incident  dose  and  r  is  the  distance 
from  the  surface  area  element  to  the  detector. 
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The  doubly  reflected  component  is  estimated  by 


= 


/ 


over  all 

illuminated  area 
area 


2  2 
'■l''2 


(5.44) 


Complexities  of  the  problem  will  almost  always 
require  the  use  of  a  computer  program  for  numeri¬ 
cal  evaluation  of  these  equations.  Cain^^  has 
shown  that  for  the  case  of  low-energy  neutrons, 
numerical  evaluation  of  these  integrals  becomes 
impractical  because  the  number  of  reflection  events 
which  must  be  considered  in  order  to  adequately 
calculate  the  penetrating  radiation  is  so  large. 
As  an  alternative,  a  method  was  developed  that 
traces  particles  through  successive  reflection 
events  by  a  Monte  Carlo  technique  which  uses  a 
random-walk  process  to  avoid  excessive  run-time 
penalties  associated  with  the  usual  Monte  Carlo 
calculation.  This  method  is  discussed  later  in 
this  section. 

In  most  albedo  calculations  of  the  transmission 
of  radiation  through  rectangular  ducts  with  bends, 
the  corner  lip  (the  corner  formed  by  the  intersection 
of  two  legs)  is  assumed  to  be  opaque  to  radiation 
in  the  main  parts  of  the  calculation  and  a  correc¬ 
tion  is  made  for  its  effect  with  a  separate  approxi¬ 
mate  calculation.  This  correction  is  usually  broken 
down  into  two  contributions:  the  corner-lip  trans¬ 
mission  effect  and  the  corner-lip  inscattering  effect. 
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Fig.  5.15.  Geometry  for  Albedo  Calculation  of  Doubly 
Reflected  Radiation  in  a  Simple  Duct  (Broad-Beam 
Source). 


The  corner-lip  transmission  effect  is  illustrated 
in  Fig.  5.16  by  the  dashed  line  that  originates  at 
the  duct  entrance  and  passes  through  the  inside 
comer.  The  radiation  penetrating  the  lip  scatters 
from  the  opposite  wall  and  contributes  to  the  dose 
at  the  detector.  In  accounting  for  this  radiation, 
the  albedo  integrations  of  Eqs.  5.43  and  5.44  can 
be  performed,  with  the  incident  dose  represented 
by  a  variable  which  depends  on  the  path  length 
in  the  duct  wall  material  of  a  ray  that  is  parallel 
to  the  initial  direction  of  the  radiation  and  termi¬ 
nates  at  the  scattering  point.  An  approximation 
may  be  made  by  assuming  that  all  particles  which 
penetrate  the  material  with  a  path  length  less 
than  a  relaxation  length  are  unattenuated  and  that 
all  particles  which  penetrate  more  than  a  relaxation 
length  are  completely  absorbed.  (LeDoux  and 
Chilton  used  the  energy  absorption  coefficient  to 
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DUCT  ENTRANCE 


Fig,  5*16#  Geometry  for  Determining  Corner-Lip  Trans¬ 
mission  Effect  in  Two-Legged  Rectangular  Duct.  This 
geometry  is  also  used  in  the  illustration  of  the  LeDoux- 
Chilton  technique  discussed  later  in  this  section. 
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Fig.  5.17.  Geometry  for  Lost-Collision  Calculation 
tangular  Duct. 


estimate  this  effective  relaxation  length  for  gamma 
rays.)  On  this  basis,  the  effect  of  the  comer-lip 
penetration  can  be  approximated  by  an  increase 
in  the  number  of  scattering  areas  beyond  the  pri¬ 
mary  scattering  areas.  The  areas  designated  as 
A^,  A^,  A^,  A^,  A^^,  and  A^^  in  Fig.  5.16  show 
some  of  these  additional  scattering  surfaces.* 

The  comer-lip  inscattering  effect  results  in  some 
of  the  radiation  being  redirected  toward  the  detec¬ 
tor  by  scattering  in  the  corner  lip.  This  component 
may  be  evaluated  by  using  a  last-collision*®  tech¬ 
nique,  for  which  the  geometry  is  shown  in  Fig. 
5.17.  The  flux  at  a  distance  in  the  corner  lip 
is  given  by 


exp 


A(E,) 


(5.45) 


*In  the  ADONIS  Monte  Carlo  calculations  reported 
in  Section  5.1.3,  Gardner  and  Mettlerl®  included  a  study 
of  the  comer-lip  effect  by  performing  two  calculations 
which  utilized  different  importance  functions.  By  proper 
selection  of  the  weighting  functions  it  was  possible 
to  eliminate  most  of  the  corner-lip  penetration  or  to 
allow  for  the  penetration.  It  was  shown  that  the  comer- 
lip  effect  can  be  important  for  some  conditions. 


of  Corner-Lip  Inscottering  Effect  in  Two-Legged  Rec- 

where  A(£p)  is  the  effective  relaxation  length  for 
the  incident  radiation  in  the  material.  The  radi¬ 
ation  per  unit  solid  angle  scattering  toward  P  by 
the  j’th  element  in  the  material  is  given  by  $  K.{E,9) 
N.  dv,  where  K.{E,9)  is  the  differential  angular 
scattering  cross  section  of  material  i  for  radiation 
incident  at  energy  E,  and  N.  is  the  atomic  density 
of  element  i  in  the  duct  wall  material.  The  flux 
at  P  due  to  radiation  of  energy  E  scattering  in  the 
lip  is  then 

<D=  /  {cfv<D„  Ie/E„,0) 

over  ^  i 

sc  attering 
volume 


X  N. 
1 


|^-[^2/A(Ei)]| 


(5.46) 


where  E^  is  the  energy  after  scattering.  This 
technique  is  analogous  to  integrating  a  point  kernel 
over  the  scattering  volume,  with  the  strength  of 
the  source  point  used  as  the  differential  scattering 
density  at  that  point. 
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Simon-Clifford  Technique  for  Cylindrical  Ducts.  — 

The  albedo  method  of  Simon  and  Clifford^®  is  one 
of  the  least  complicated  of  the  methods  available 
for  calculating  the  transmission  of  thermal  and 
intermediate-energy  neutrons  through  cylindrical 
ducts.  It  considers  only  the  neutrons  entering  the 
open  end  of  the  duct,  and  transmission  is  assumed 
to  be  either  by  line  of  sight  or  by  scattering  from 
the  duct  walls.  The  wall  scattering  is  calculated 
by  using  a  single-energy  “spectrum-averaged”  ma¬ 
terial  albedo. 

In  the  case  of  a  straight  cylindrical  duct,  the 
albedo,  A^,  is  defined  as  the  fraction  of  incident 
neutrons  which  are  reradiated  (see  Chapter  4).  In 
this  form  of  the  albedo,  the  reflected  radiation 
has  been  integrated  over  all  exit  directions,  the 
reflected  angular  distribution  being  expressed  in 
the  general  form  ^  +  (2y  cos  d)/2TT,  where  is 
the  fraction  reradiated  isotropically,  y  is  the  frac¬ 
tion  reradiated  with  a  cosine  distribution,  and 
/3  +  y  =  1.  The  flux  at  distance  Z  along  a  duct 
of  radius  a  is  then  given  by 
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Fig.  5.18.  Geometry  for  Calculating  Neutron  Trans¬ 
mission  Through  Cylindrical  Duct  by  Simon-Clifford 
Method. 

where  A'^  is  an  empirically  determined  albedo  pa¬ 
rameter. 

If  a  series  of  bends  are  involved,  the  flux  at  a 
distance  Z  from  the  last  bend  may  be  expressed  as 


(5.47) 


where  ^Z)  is  in  the  same  units  as  the  surface 
source  strength 

The  first  term  in  this  equation  is  the  line-of-sight 
component  given  in  Section  5.1.1.  The  second 
term,  which  accounts  for  the  wall-scattering  com¬ 
ponent,  is  a  strong  function  of  the  material  albedo. 
Since  A^  is  of  the  order  of  0.1  for  fast  neutrons 
and  the  second  term  is  small  for  long  thin  ducts 
(a  «  Z),  Eq.  5.47  simplifies  to  that  for  the  uncol¬ 
lided  flux.  However,  for  thermal  neutrons,  which 
have  an  albedo  of  about  0.8,  the  second  term 
becomes  important,  and  the  scattered  flux  can 
dominate. 

This  approach^®  was  extended  to  cover  the  case 
of  a  bent  cylindrical  duct  (see  Fig.  5.18)  by  assum¬ 
ing  that  the  albedo  (in  the  area  of  the  bend)  of 
the  neutrons  transmitted  through  the  first  leg  of 
the  duct  is  the  source  for  calculation  of  transmis¬ 
sion  through  the  second  leg.  When  both  legs  are 
of  equal  diameter  and  their  angular  separation  is 
6,  the  total  transmission  at  a  distance  Z  ^  along  the 
second  leg  is  given  by 


/S  +  2y  sin  d  ^ 
sin  6, 


/3  +  2y  sin 


(5.49) 


where  the  duct  segments  are  of  lengths  L^,  L^,  .  . . 
and  are  separated  by  angles  6^,  9^,  ...  0„_j. 

Ducts  which  follow  a  smooth  curve  through  the 
shield  may  be  treated  by  the  above  method  pro¬ 
vided  that  the  radius  of  curvature  is  large  with 
respect  to  the  radius  of  the  duct.  ^  The  curved 
path  may  be  divided  into  a  series  of  equivalent 
straight  sections  of  length  equal  to  the  maximum 
chord  which  can  be  drawn  internal  to  the  duct. 
The  angular  separation  between  successive  chords 
may  then  be  used  as  the  angle  between  the  equiv¬ 
alent  sections. 

Horton*®  treated  the  case  of  helical  ducts  with 
an  albedo  approach  similar  to  the  Simon-Clifford 
method.  For  a  helical  path  divided  into  n  equiv¬ 
alent  straight  sections  of  radius  a  and  length  / 
deviating  successively  by  constant  angle  ip,  the 
flux  is  approximated  by 


0(Z,) 


/3  +  2y  sin  9 
sin  9 


(5.48) 
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This  may  be  compared  with  Eq.  5.49  which  would 
be  the  equivalent  Simon-Clifford  expression  provided 
that  all  values  of  6  equal  ifi  and  all  values  of  L 
equal  /. 

The  usefulness  of  the  semiempirical  albedo  tech¬ 
niques  described  above  depends  on  the  availability 
of  suitable  values  of  the  albedos  and  angular  dis¬ 
tributions  of  the  reflected  neutrons. 

LeDoux-Chilton  Technique  for  Rectangular  Ducts. 
A  method  of  calculating  gamma-ray  transmission 
through  two-legged  rectangular  ducts  was  formu¬ 
lated  by  LeDoux  and  Chilton,^®  and  the  same 
general  technique  was  used  by  Song*®  for  neutron 
transmission.  Equations  were  derived  for  calculat¬ 
ing  radiation  stream.ing  down  two-legged  concrete 
ducts  using  as  transport  mechanisms;  (1)  a  single¬ 
reflection  albedo  for  surfaces  that  can  be  “seen” 
by  both  the  source  and  the  detector,  and  (2)  scat¬ 
tering  by  material  in  the  comer  lip  of  the  duct. 

In  the  absence  of  sufficient  differential  albedo 
data  at  the  time  of  this  work,  a  total  albedo^** 
was  used  and  the  radiation  was  assumed  to  emerge 
isotropically.  To  simplify  the  equations,  it  was 
assumed  that  all  points  on  a  major  scattering  area 
(numbered  areas  of  Fig.  5.16)  could  be  considered 
to  have  the  same  angle  of  incidence  from  the 
source  and  the  same  reflection  angle  to  the  detec¬ 
tor.  Formulas  were  then  derived  for  the  detector 
response  due  to  the  first-order  scattering  from 
each  of  the  surfaces.  Later  work^*~^®  with  more 
accurate  differential  albedo  data  showed  that  this 
method  tends  to  underpredict  because  multiple 
reflections  are  neglected.  That  the  method  pre¬ 
dicts  experimental  data  as  well  as  it  does  is 
attributed  to  a  compensating  overprediction  caused 
by  the  assumption  of  isotropic  scattering. 

Chapman^'*  extended  the  LeDoux-Chilton  approach 
to  include  double  reflections  and  second-order 
effects  such  as  a  wall  backscatter  followed  by  a 
corner-lip  inscatter.  In  place  of  the  total  albedos 
with  the  assumption  of  isotropy  of  the  reflected 
radiation  used  by  LeDoux  and  Chilton,  Chapman 
substituted  the  semiempirical  formula  for  the  differ¬ 
ential  dose  albedo  that  was  derived  by  Chilton 
and  Huddleston^®  (see  Section  5.1.5).  The  two 
parameters  in  this  formula  were  determined  by  a 
least-squares  fit  to  the  Monte  Carlo  data  of  Raso.^® 

The  complexity  of  the  interaction  combinations 
considered  (e.g.,  one  backscatter,  one  inscatter, 
one  backscatter  plus  one  penetration,  two  back- 
scatters,  and  one  backscatter  plus  one  inscatter) 


dictated  that  Chapman  use  four  computer  programs 
to  calculate  the  following: 

1.  the  original  LeDoux-Chilton  formulas  (singly 
reflected  component), 

2.  the  LeDoux-Chilton  formulas  with  the  scattering 
areas  broken  into  smaller  increments, 

3.  the  contribution  due  to  second-order  scattering 
from  the  duct  walls, 

4.  the  contribution  due  to  inscatter  by  the  corner 
lip  of  the  duct  of  either  direct  or  first-order 
wall-scattered  radiation. 

The  singly  reflected  component  was  evaluated  by 
integrating  over  each  of  the  numbered  areas  in  Fig. 
5.16.  In  evaluating  penetration  through  the  corner 
lip,  the  assumption  was  made  that  rays  penetrating 
less  than  one  relaxation  length  of  material  are 
unattenuated  and  that  those  encountering  more  than 
one  relaxation  length  of  material  do  not  penetrate. 
Scattering  in  the  corner  lip  was  calculated  using 
the  singly  reflected  formula  which  LeDoux  and 
Chilton  derived  on  the  basis  of  Klein-Nishina  dif¬ 
ferential  scattering  probabilities. 

To  perform  the  double-reflection  calculations. 
Chapman  divided  the  duct  into  a  larger  mesh  con¬ 
sisting  of  12  scattering  areas  as  shown  in  Fig. 
5.19,  and  combinations  which  would  allow  transport 
from  the  source  to  the  detector  with  two  scatters 
were  considered.  Combinations  whose  percentage 
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Fig,  5,19.  Scattering  Areas  for  Albedo  Calculation  of 
Doubly  Reflected  Gamma  Rays  in  Two-Legged  Rectangu¬ 
lar  Duct. 
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contribution  to  the  total  dose  was  considered  to 
be  small  were  deleted.  Specifically  these  were 
reflection  from  surface  1  to  surfaces  10  and  12 
and  reflection  from  surfaces  2  and  4  to  surfaces 
10,  11,  and  12.  The  incidence  of  singly  reflected 
radiation  on  the  corner  lip  of  the  duct  was  com¬ 
bined  with  the  lip  inscatter  calculation  to  deter¬ 
mine  the  contribution  due  to  one  backscatter  plus 
one  inscatter.  In  all  cases  the  energy  degradation 
in  a  reflection  was  assumed  to  be  that  of  a  single 
Compton  scattering. 

In  calculating  neutron  transmission  through  rec¬ 
tangular  ducts,  Song'^  used  Chapman’s  modifi¬ 
cation  of  the  LeDoux-Chilton  method  with  the  ex¬ 
ception  that  the  reflected  neutron  spectrum  was 
assumed  to  consist  of  two  neutron  energy  groups, 
the  number  and  energy  in  each  group  being  based 
on  the  Monte  Carlo  data  of  Allen,  Futterer,  and 
Wright.^®  Song  derived  a  semiempirical  formula 
for  use  in  the  calculations  in  terms  of  a  single¬ 
energy  parameter  for  the  differential  neutron  dose 
albedo  in  a  manner  analogous  to  that  used  by 
Chilton  and  Huddleston^®  in  their  gamma-ray  cal¬ 
culations.  Values  for  the  energy-dependent  param¬ 
eter  were  obtained  by  a  least-squares  analysis 
that  gave  the  best  fit  to  the  data  of  Allen  et  al. 

Monte  Carlo  Methods  for  Rectangular  Ducts.  — 
Generally,  experiments  and  calculations  have  shown 
that  gamma  rays  produced  by  neutron  absorptions 
in  a  concrete  shield  surrounding  a  rectangular 
duct  make  a  significant  contribution  to  the  total 
dose  in  the  duct,  and  analysis  of  this  component 
requires  a  knowledge  of  the  distribution  of  low- 
energy  neutrons  throughout  both  the  duct  and  the 
shield.  These  distributions  can  be  determined  ex¬ 
actly  by  analog  Monte  Carlo  calculations,  but,  as 
was  pointed  out  in  Section  5.1.3,  the  machine 
time  required  to  trace  neutrons  until  they  approach 
thermal  equilibrium  can  become  excessive.  To 
circumvent  this,  Cain  developed  a  technique*^  that 
uses  a  random-walk  Monte  Carlo  approach  in  which 
particles  are  traced  by  selecting  random  path  direc¬ 
tions  inside  the  duct  while  prohibiting  particle 
penetration  into  the  shielding  material.  Thus, 
tracking  through  the  large  number  of  interactions 
prior  to  emergence  from  the  wall  material  is  averted 
and  no  time  is  wasted  in  tracking  particles  which 
never  emerge.  Each  time  a  particle  encounters 
the  wall,  it  is  assumed  to  re-emerge  with  a  reduced 
weight,  given  by  the  incident  weight  times  the 
albedo  of  the  wall.  Fluxes  at  various  points 
along  the  duct  are  then  calculated  either  by  mak¬ 


ing  a  statistical  estimate  from  each  encounter  to 
point  detectors  or  by  determining  the  track  length 
per  unit  volume  in  a  detector  region.  The  most 
satisfactory  detector  was  found  to  be  a  thin  region 
extending  across  the  duct  and  perpendicular  to  the 
duct  axis. 

The  method  was  modified  by  Maerker  and  Cain^^ 
to  include  doubly  differential  thermal-neutron  cur¬ 
rent  and  fast-neutron  dose  albedo  data  (differential 
in  reflected  solid  angle  and  energy)  previously 
determined^®"®”  by  Monte  Carlo  calculations  for 
infinite  slabs  of  concrete  (see  Chapter  4).  In 
addition,  a  provision  was  included  to  calculate  the 
capture-gamma-ray  dose  rate  arising  from  neutron 
capture  in  the  duct  walls. 

In  the  modified  code  the  incident  neutrons  are 
reflected  with  reduced  “weight”  into  various  energy 
groups  and  directions,  the  probabilities  for  which 
are  predetermined  from  the  doubly  differential 
albedo  data.  The  neutrons  are  followed  in  the 
random-walk  procedure  until  they  either  escape  out 
the  front  or  rear  end  of  the  duct  or  are  killed  by 
Russian  roulette  when  their  weight  falls  below  a 
predetermined  value.  Statistical  estimations  of  the 
fluxes  to  point  detectors  are  made  to  obtain  the 
results. 

The  best  accuracy  is  obtained  when  the  calcu¬ 
lations  are  done  separately  for  three  energy  ranges: 
(1)  Neutrons  with  energies  from  200  keV  to  8  MeV 
are  divided  into  six  energy  groups  and  are  sampled 
to  obtain  the  fast  flux  contribution  to  the  dose 
rate.  (2)  The  dose  rate  from  neutrons  of  inter¬ 
mediate  energy  (0.5  eV  to  200  keV),  including  all 
neutrons  that  have  slowed  down  to  intermediate 
and  thermal  energies  from  higher  energy  groups, 
is  determined  with  a  13-group  calculation  that 
involves  sampling  the  entire  energy  range  from 
thermal  to  8  MeV  but  scoring  only  those  neutrons 
with  energies  below  200  keV.  (3)  The  dose  rate 
from  source  thermal  neutrons  (subcadmium,  <0.5 
eV)  is  then  determined  in  a  separate  calculation 
using  a  single-velocity  model. 

Different  Monte  Carlo  albedo  data  are  used  in 
each  energy  range.  The  albedos  used  for  the  fast- 
neutron  calculation  are  those  reported  for  concrete 
by  Maerker  and  Muckenthaler;  they  do  not  cover 
neutrons  reflected  with  energies  less  than  200 
keV.  The  albedos  used  for  the  intermediate-energy 
neutrons  are  those  determined  by  Coleman  et  al.^^ 
The  latter  include  the  probability  of  higher  energy 
neutrons  being  reflected  as  thermal  neutrons.  The 
albedos  for  neutrons  that  have  thermal  energy  both 
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before  and  following  reflection  are  those  deter¬ 
mined  by  Maerker  and  Muckenthaler^°  using  a 
one-velocity  model.  Finally,  the  equivalent  of  an 
albedo  is  used  to  calculate  the  capture-gamma- 
ray  contribution  to  the  dose,  the  results  being 
correlated  on  the  basis  of  the  gamma-ray  dose 
emitted  per  steradian  per  incident  neutron  at  a 
particular  angle  of  incidence  and  within  a  certain 
energy  group.  The  capture-gamma-ray  albedos  for 
incident  thermal  neutrons  were  developed  by  Maerker 
and  Muckenthaler^®  and  those  for  incident  inter¬ 
mediate  neutrons  by  Coleman  et  at. 

In  this  modified  method  no  neutron  is  allowed  to 
penetrate  a  corner  wall  during  the  course  of  its 
random  walk,  although  a  statistical  estimate  of  the 
flux  that  includes  a  corner-lip  penetration  in  the 
last  flight  may  be  made.  This  is  accomplished 
simply  by  neglecting  the  contribution  when  the 
path  through  a  corner  is  greater  than  1  mean  free 
path  and  including  the  unattenuated  contribution 
if  the  path  is  less  than  1  mean  free  path.  No  con¬ 
sideration  is  given  to  the  comer-lip  inscattering 
effect. 

Comparison  with  Experiment.  —  The  accuracy  of 
albedo-type  duct  streaming  calculations  depends 
on  the  accuracy  of  the  albedo  data  used,  on  the 
detail  with  which  the  scattering  surface  mesh  and 
the  energy  groups  are  defined,  and  on  the  selection 
of  the  multiply  scattered  components  to  be  included. 
The  comparisons  given  in  the  following  paragraphs 
show  the  extent  of  agreement  that  has  been  ob¬ 
tained  when  the  albedo  approach  has  been  used 
to  analyze  a  variety  of  experiments.  The  reason 
that  most  of  the  comparisons  are  for  gamma  rays 
from  isotope  sources  attenuated  through  bent  rec¬ 
tangular  ducts  is  that  a  number  of  experiments 
have  been  performed  with  such  sources  in  con¬ 
junction  with  studies  on  the  design  and  evaluation 
of  shelters  to  protect  against  fallout  radiation. 
Only  typical  data  have  been  selected  for  the  com¬ 
parisons  presented  here,  since  an  exhaustive  sum¬ 
mary  has  been  published  by  Huddleston  and  Wilcox- 
son.  ^ ' 

Comparisons  with  51  different  measurements  were 
made  by  Chapman,  who  used  the  extended  LeDoux- 
Chilton  technique  to  calculate  gamma-ray  dose 
rates  at  the  exit  end  of  two-legged  L-shaped  rec¬ 
tangular  ducts,  the  source  in  each  case  being  a 
gamma-ray-emitting  isotope  at  the  mouth  of  the  duct. 
The  results  are  presented  in  Table  5.4,  where  W, 
Lj,  and  are  as  defined  in  Fig.  5.16.  Agreement 


is  generally  within  30%  for  a  ®°Co  or  ^‘‘Na  source, 
within  45%  for  ^^^Cs,  and  within  a  factor  of  2 
for  ^®®Au.  The  analysis  overpredicted  the  ex¬ 
perimental  results  fairly  consistently,  but  not  by 
a  consistent  margin. 

Chapman  also  calculated  the  gamma-ray  dose  rates 
along  the  axes  of  two-  and  three-legged  rectangular 
ducts  and  some  of  them  are  compared  in  Figs. 
5.20  through  5.23  with  the  dose  rates  measured 
at  several  points  along  the  duct  axes  for  three 
sources  of  different  energy. Figures  5.21 
and  5.22  also  show  comparisons  with  a  simplified 
formula  (Eq.  5.51)  proposed  by  Ingold  and  Huddles¬ 
ton,^®  which  is  discussed  in  Section  5.1.5. 

As  is  apparent  from  Table  5.4,  a  number  of  the 
experiments  with  which  Chapman  made  his  com¬ 
parisons  were  performed  by  Terrell  et 
Although  not  considered  by  Chapman  in  his  com¬ 
parisons,  some  of  the  experiments  performed  by 
Terrell  et  al.^^  were  to  investigate  the  effect  of 
geometry  and  source  energy  on  the  attenuation  of 
gamma  rays  through  rectangular  concrete  ducts 
with  two  right-angle  bends,  the  results  of  which 
show  trends  which  can  be  anticipated  in  calcu- 
tions.  They  found  that  for  6-ft-square  concrete 
ducts  it  made  little  difference  whether  the  gamma- 
ray  source  was  ®°Co  or  ^®^Cs  or  whether  the  right- 
angle  bends  formed  a  U  shape  (Fig.  5.24a)  or  a  Z 
shape  (Fig.  5.246).  These  data,  plotted  in  Fig. 
5.24c,  show  a  maximum  deviation  of  a  factor  of 
1.5  when  properly  normalized.  The  effect  of  geom¬ 
etry  and  energy  becomes  more  noticeable  in  the 
case  of  a  1-ft-square  duct.  Typical  examples  of 
the  effect  of  gamma  source  energy  and  duct  shape 
are  shown  in  Figs.  5.25  and  5.26  respectively. 

In  his  studies  Chapman  also  calculated  the  con¬ 
tribution  of  a  given  scattering  area  on  the  gamma- 
ray  dose  rate  in  an  L-shaped  3-ft-square  concrete 
duct.  The  scattering  area  chosen  was  the  one 
shown  as  Sj  in  Fig.  5.27a,  it  being  assumed  that 
gamma  rays  from  a  ®°Co  source  at  the  mouth  of  the 
duct  were  scattered  from  this  area  to  the  detector 
at  the  exit  end  of  the  duct.  The  dose  rate  due  to 
gamma  rays  initially  scattered  from  was  calcu¬ 
lated  to  be  0.379  mr/hr,  which  compares  quite 
favorably  with  a  measured  dose  rate  of  0.404  mr/hr. 

During  the  experiment.  Chapman  measured  the 
gamma-ray  spectrum  at  the  exit  end  of  the  duct  both 
with  and  without  the  scattering  area  being  shad¬ 
owed  from  the  source.  The  results  are  shown  in 
Fig.  5.276.  Shielding  of  had  a  more  pronounced 
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Table  5.4.  Comparison  of  Gamma-Ray  Dose  Rates  Predicted  by  LeDoux-Cbilton  Technique 
with  Dose  Rates  Measured  at  Exit  End  of  Two-Legged  Duct® 


Gamma-Ray 

Source 

w 

L^/W 

Dose  Rate  (mr/hr) 

%  Difference*’ 

Reference 

Calculated 

Measured 

0.34-curie 

11  in. 

1.90 

1.65 

87.3 

125 

-30 

Green^ 

2.06 

44.5 

61 

-27 

2.46 

27.1 

30.5 

-11 

3.68 

8.46 

7.31 

-tie 

3.58 

2.06 

6.17 

7.3 

-15 

2.86 

2.61 

2.7 

-3 

3.68 

1.30 

1.3 

0 

0.6-curie 

11.1  in. 

3.54 

1.73 

17.4 

15.6 

-til 

Eisenhauer"* 

2.79 

4.94 

3.7 

-t-33 

3.51 

2.65 

2.02 

-131 

55-curie 

12  in. 

3.50 

2.0 

916 

852 

-18 

Terrell® 

3.0 

317 

243 

-130 

4.0 

140 

110 

+28 

2.4-curie  ®°Co 

3  ft 

2.0 

1.67 

20.6 

17.5 

+18 

Chapman* 

2.0 

12.6 

12.1 

+4 

2.34 

8.35 

7.1 

+18 

2.5 

1.50 

14.5 

13.5 

+7 

1.83 

8.42 

9.1 

-8 

2.0 

6.70 

6.4 

+5 

2.5 

3.79 

3.7 

+2 

3.67-curie 

6  ft 

1.33 

1.83 

15.4 

11.8 

+31 

Terrell® 

2.50 

6.56 

4.75 

+38 

3.17 

3.47 

2.42 

+43 

1.66 

1.83 

7.85 

7.30 

+8 

2.50 

3.46 

2.73 

+27 

3.17 

1.85 

1.39 

+33 

2.0 

1.83 

4.71 

4.56 

+3 

2.50 

2.12 

1.79 

+  18 

3.17 

1.14 

0.935 

+21 

80-curie 

12  in. 

3.5 

2.0 

606 

430 

+41 

Terrell® 

3.0 

208 

132 

+58 

4.0 

90 

90 

+41 

6  ft 

2.17 

1.83 

36.5 

35.5 

+3 

2.33 

19.7 

19.6 

0 

1.52-curie 

6  ft 

2.0 

1.83 

0.858 

0.714 

+20 

3.17 

0.207 

0.186 

+11 

15.5-curie  ^^^Au 

3  ft 

2.0 

1.67 

55.8 

37.6 

+48 

Chapman^ 

2.0 

3.4 

23.3 

+46 

2.67 

15.7 

11.1 

+41 

11  in. 

4.125 

2.04 

85 

46.7 

+83 

2.86 

35.2 

19.0 

+85 

3.27 

24 

13.5 

+78 

25 
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Table  5.4  (continued) 


Gamma  "Ray 
Source 


Dose  Rate  (mr/hr) 

W  L  /W  L  /W  - 

^  Calculated  Measured 


%  Difference*’ 


Reference 


1.66 

1.83 

8.77* 

6.78 

4-29 

2.50 

3.84 

2.80 

+37 

3.17 

2.05 

1.50 

+37 

2.17 

1.83 

4.17 

3.64 

+15 

2.50 

1.88 

1.67 

+13 

3.17 

1.02 

0.912 

+11 

2.84 

1.83 

2.02 

1.94 

+9 

2.50 

0.931 

0.828 

+12 

3.17 

0.475 

0.462 

+3 

Terrell^ 


®Table  taken  from:  J.  M.  Chapman,  Computer  Calculations  of  Dose  Rates  in  Two-Legged  Ducts  Using  the  Albedo 
Concept,  Naval  Civil  Engineering  Laboratory  Report  NCEL-R-264  (Oct.  24,  1963). 

^Percent  difference  =  (calculated  —  measured/measured)  100. 

‘’D.  W.  Green,  Attenuation  of  Gamma  Radiation  in  a  Two-Legged  11-inch  Rectangular  Duct,  Naval  Civil  Engi¬ 
neering  Laboratory  Report  NCEL-R-195  (May  2,  1962). 

‘*C.  Eisenhauer,  Scattering  of  Cobalt-60  Gamma  Radiation  in  Air  Ducts,  National  Bureau  of  Standards  Technical 
Note  74  (1960). 

®C.  W.  Terrell,  A.  J.  Jerri,  and  R.  O.  Lyday,  Jr.,  Radiation  Streaming  in  Ducts  and  Shelter  Entranceways,  Armour 
Research  Foundation  Report  ARF  1158-A02-7  (April  1962). 

*J.  M.  Chapman,  Gamma  Dose  Rates  and  Energy  Spectra  in  a  3-Foot  Square  Duct,  Naval  Civil  Engineering  Lab¬ 
oratory  Report  NCEL-N-443  (June  1962). 

^C.  W.  Terrell  et  al.,  Radiation  Streaming  in  Shelter  Entranceways ,  Armour  Research  Foundation  Report  ARF 
1158-12  (October  1960). 

*’J.  M.  Chapman,  The  Variation  of  Dose  Attenuation  of  Two-Legged  Concrete  Ducts  with  Incident  Gamma-Ray 
Energy,  Naval  Civil  Engineering  Laboratory  Report  NCEL-TN-707  (April  1965). 

*The  values  for  ^'^Na  are  the  sum  of  the  values  obtained  using  initial  gamma-ray  energies  of  1.37  and  2.75  MeV. 

>0.  W.  Terrell  and  A.  J.  Jerri,  Radiation  Streaming  in  Shelter  Entranceways,  Armour  Research  Foundation  Report 
ARF  1158-A01-5  (July  1961). 


effect  on  the  higher  energy  radiation  than  on  the 
multiply  reflected  low-energy  component.  This 
effect  would  be  expected  since  the  percentage  re¬ 
duction  in  the  single-reflected  area  that  results 
from  blocking  Sj  is  greater  than  the  percentage 
reduction  in  possible  combinations  of  multiple- 
reflection  areas. 

In  a  subsequent  experiment.  Chapman  and  Grant 
investigated  the  relative  effectiveness  of  copla- 
nar  and  noncoplanar  three-legged  rectangular  ducts. 
For  the  cases  investigated  (11-in. -square  duct, 
®°Co  source,  =  30  in.,  =  31.5  in.,  =  36 
in.),  the  attenuation  provided  by  the  noncoplanar 
configuration  was  found  to  be  greater  by  a  factor 
of  2  than  that  provided  by  the  coplanar  case. 
This  effect  cannot  be  explained  by  analysis  be¬ 
cause  third  or  higher  order  reflections  are  diffi¬ 
cult  to  treat  with  the  albedo  theory.  The  differ¬ 


ence  could  be  due  to  geometric  effects  on  the 
scattering  angles. 

Clifford^’  used  differential  albedos to  cal¬ 
culate  dose  rates  corresponding  to  those  measured 
in  experiments  designed  to  determine  the  impor¬ 
tance  of  multiple  reflection  on  the  transmission 
of  *®^Cs  gamma  rays  through  8-in.-,  12-in.-,  and 
3-ft-square  concrete  ducts.  In  the  experiments, 
collimated  gamma  rays  were  incident  on  the  duct 
entrance  at  an  angle  to  the  duct  axis  such  that 
the  area  in  which  the  first  reflection  could  occur 
was  near  the  entrance  of  the  duct  (see  Fig.  5.28a). 
The  dose  rate  along  the  duct  center  lines  was 
first  measured  with  the  concrete  duct  unmodified 
(Condition  j4)  and  then,  in  order  to  estimate  the 
influence  of  each  reflection  area  in  the  duct,  with 
a  lead  lining  covering  all  inside  surfaces  of  the 
duct  except  the  first-scatter  area  (Condition  B) 
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Fig.  5.20.  Comparison  of  Albedo  Calculations  with 
Measurements  of  Gamma-Ray  Dose  Rates  Along  Axis  of 
l-ft-square  U-Shaped  Concrete  Duct:  80-curie  ^^^Cs 
Gamma-Ray  Source.  (Experimental  data  from  Terrell 
et  al,,  ref.  32;  calculations  from  Chapman,  ref.  24.) 
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Fig.  5.22.  Comparison  of  Albedo  Calculations  with 
Measurements  of  Gamma-Ray  Dose  Rates  Along  Axis 
of  6-ft-square  L-Shaped  Concrete  Duct:  3.67-curie 

^^Co  Source.  The  results  designated  ‘‘simplified  for¬ 
mula**  correspond  to  the  Ingold-Hudd  leston  equation 
given  In  Section  5.1.5,  (Experimental  data  from  Terrel! 
et  al.f  ref.  33;  calculations  from  Chapman,  ref,  24.) 
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Fig*  5.21.  Comparison  of  Albedo  Calculations  with 
Measurements  of  Gamma-Ray  Dose  Rates  Along  Axis  of 
6-ft-square  L-Shaped  Duct:  3,67-curie  ^^Co  Source* 
The  results  designated  “simplified  formula**  correspond 
to  the  Ingold-Hudd leston  equation  given  in  Section  5.1 .5, 
(Experimento I  dota  from  Terrell  et  a/.,  ref.  33;  colculo- 
tions  from  Chopman,  ref,  24.) 
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Fig.  5.23.  Comparison  of  Albedo  Calculations  with 
Measurements  of  Gamma-Ray  Dose  Rates  Along  Axis  of 
6-ft-square  L-Shaped  Duct:  4.2-curie  ^^Na  Source. 

(Experimental  data  from  Terrell  ef  al.,  ref,  34;  calcula¬ 
tions  from  Chapman,  ref.  24-) 
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Fig.  5.28a.  Schematic  of  Square  Concrete  Duct  Show¬ 
ing  First  Scatter  Area. 
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Fig.  5,28b.  Measurements  and  Single-Reflection  Anal¬ 
ysis  of  Gamma-Ray  Dose  Rates  Along  Axis  of  8-in.- 
square  Concrete  Duct  With  and  Without  Lead  Linings: 
^^^Cs  Source.  {F  rom  Clifford,  ref.  21.) 


and  also  with  a  lead  lining  in  the  entire  duct 
(Condition  C).  Figure  5.286  shows  the  results 
of  the  measurements  for  all  three  conditions  in 
the  8-in.  duct,  together  with  a  single-reflection 
analysis  based  on  the  differential  albedo  for  Con¬ 
dition  B,  and  Fig.  5.29  shows  the  fraction  of  the 
total  dose  which  is  due  to  reflection  from  sur¬ 
faces  not  exposed  to  the  primary  radiation  for 
all  the  ducts.  It  was  apparent  from  the  results 
of  all  three  duct  sizes  that  the  multiply  reflected 
component  increases  to  become  as  much  as  20 
to  30%  of  the  total  dose  as  the  duct  length  is 
increased. 

Brodeur  and  Batter^®  used  the  transmission 
and  backscatter  data  of  Raso^^  to  predict  the 
effect  of  a  protruding  air  vent  on  the  dose  rate 
measured  in  an  underground  shelter  that  was  due 
to  a  ®°Co  source  above  the  shelter.  The  pre¬ 
dictions  were  made  to  compare  with  measurements 
made  for  the  geometry  shown  in  Fig.  5.30a.  The 
experimental  results  are  shown  in  Fig.  5.306, 
where  they  are  plotted  as  the  dose  rate  at  the 
detector  position  with  the  protruding  cylinder 
(D  ,  )  minus  the  dose  rate  without  the  cylinder 
(D  /  )  divided  by  the  dose  rate  incident  on  the 
cylinder  (Dp.  The  detector  was  located  in  sev¬ 
eral  positions  below  the  air  vent  and  various 
cylinder  heights  were  used.  The  calculations 
predicted  dose  rates  that  were  lower  by  factors 
of  2  to  5  than  those  observed  experimentally. 
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Fig.  5.29.  Ratio  of  Multiply  Scattered  Gamma-Ray 
Dose  Rate  to  Total  Dose  Rate  Along  Axes  of  8-in.-, 
12-in.-,  and  3-ft-square  Concrete  Ducts;  '^^Cs  Source. 
(From  Clifford,  ref.  21.) 

It  was  pointed  out  above  that  Song^^  ap¬ 
plied  the  LeDoux-Chilton  albedo  technique  to 
the  problem  of  neutron  streaming  in  two-legged 
rectangular  concrete  ducts,  the  calculational  model 
being  identical  in  geometric  detail  with  that  used 
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Fig.  5.30a.  Schematic  of  Underground  Shelter  and 
Protruding  Air  Vent. 


by  Chapman  for  gamma-ray  streaming  calcu¬ 
lations.  The  duct  was  assumed  to  be  3  ft  square 
in  cross  section,  and  the  sources  considered  were 
14.7-  and  2.5-MeV  neutrons.  The  analysis  in¬ 
cluded  all  single  and  double  scattering  events 
and  used  a  fit  to  the  Monte  Carlo  data  of  Allen 
et  al.^^  to  determine  the  energy  and  angular  de¬ 
pendence  of  the  neutron  dose  albedo.  A  com¬ 
parison  of  the  results  with  measurements  made 
by  Doty^^  (Fig.  5.31)  shows  good  agreement.* 
(The  solid  line  represents  the  calculations  for 
both  the  2.5-  and  the  14.7-MeV  sources  since 
there  was  very  little  difference  between  the  two 
sets  of  results.) 

Horton  and  Halliday^®  used  the  Simon-Clifford 
albedo  method  to  calculate  the  transmission  of 
neutrons  through  straight  cylindrical  ducts.  The 
calculations  corresponded  to  an  experiment  in 
which  1-,  2-,  and  4-in.-diam  steel-lined  ducts 
penetrated  a  water  shield  that  was  adjacent  to 
an  extended  plane  source  of  thermal  neutrons. 
A  current  albedo,  /S ,  of  0.55  was  used,  and  a 
cosine  distribution  was  assumed.  The  comparison 
shown  in  Fig.  5.32  demonstrates  that  the  Simon- 
Clifford  equation  gives  good  agreement  with  ex¬ 
periments  in  which  the  transmission  of  thermal 
neutrons  through  cylindrical  ducts  is  measured. 


*As  was  discussed  in  Section  5.1.3,  Gardner  and 
Mettler^®  used  the  ADONIS  analog  Monte  Carlo  code 
to  perform  calculations  corresponding  to  similar  meas¬ 
urements  by  Doty  12  for  a  14.7-MeV  source,  and  they 
compared  these  results  with  albedo  calculations  made 
by  Songl^  and  by  Maerker.  1“* 
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Fig.  5.30f).  Measurements  of  Dose  Rates  in  an  Under¬ 
ground  Shelter  Resulting  from  Gamma  Rays  Entering 
Through  a  Protruding  Air  Vent:  ^®Co  Gamma-Ray  Source. 
The  solid  lines  are  smooth-curve  fits  to  the  doto.  (From 
Brodeur  and  Batter,  ref.  38.) 


provided  that  reasonable  assumptions  are  made 
for  the  thermal-neutron  albedos. 

Cain  *  ^  performed  a  series  of  calculations  with 
his  original  Monte  Carlo  albedo  code  that  cor¬ 
responded  to  experiments  in  which  the  ORNL 
Tower  Shielding  Reactor  II  was  used  as  the  source 
and  the  thermal-neutron  fluxes  were  measured 
in  a  three-legged  rectangular  concrete  tunnel  con¬ 
necting  two  underground  bunkers.  The  bunkers 
were  12-ft  cubes  connected  by  a  3-ft-wide  by 
8-ft-high  tunnel  as  shown  in  Fig.  5.33a.  Two 
sets  of  measurements  were  made,  one  originating 
at  the  center  line  of  the  top  bunker,  in  which  case 
the  top  of  the  top  bunker  was  left  open,  and  the 
other  originating  at  the  center  line  of  the  front 
bunker,  in  which  case  the  front  of  the  front  bunker 
was  left  open.  The  calculations  assumed  that 
the  total  current  albedo  for  thermal  neutrons  was 
0.8  and  that  the  reflected  neutrons  had  a  cosine 
distribution.  The  comparison  given  in  Fig.  5.33b 
shows  that  for  this  configuration  good  agreement 
was  obtained  between  the  calculations  and  the 
experiment.  Cain  also  used  the  albedo  Monte 
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Carlo  code  to  calculate  thermal-neutron  trans¬ 
mission  through  cylindrical  ducts,  and  obtained 
results  that  were  in  good  agreement  with  those 
determined  by  the  Simon-Clifford  method. 

Maerker  and  Muckenthaler  carried  out  an  ex¬ 
tensive  series  of  experiments  and  calculations^®’^® 
to  test  the  accuracy  of  the  Monte  Carlo  albedo 
method  as  programmed  in  the  AMC  code.^^  Three 
configurations  of  a  3-ft-square  duct  were  used: 
a  straight  duct  45  ft  long,  a  two-legged  duct  with 
a  right-angle  bend  located  15  ft  down  the  first 
leg,  the  second  leg  being  30  ft  in  length,  and  a 
three-legged  duct  with  two  right-angle  bends,  the 
two  bends  being  located  15  ft  down  the  first  and 
second  legs  and  the  third  leg  being  20  ft  in  length. 
In  the  experiments  a  neutron  beam  from  the  Tower 
Shielding  Reactor  II  entered  each  duct  through 
an  approximately  1-ft^  area  located  at  the  geo- 


Fig.  5.31.  Comparison  of  Albedo  Calculations  with 
Measurements  of  Neutron  Dose  Rates  in  Second  Leg  of 
3-ft-square  Two-Legged  Concrete  Duct:  2,5-  and  14,7- 
MeV  Neutron  Sources,  In  the  experiment  the  source  was 
1  ft  outside  the  duct  and  =  15  ft;  in  the  calculation 
the  source  was  assumed  to  be  at  the  duct  entrance  and 
=  16  ft.  (From  Song,  ref.  13.) 


metric  center  of  the  duct  mouth.  The  neutrons 
were  incident  at  an  angle  of  45°  to  one  of  the 
duct  walls.  This  angle  of  incidence  was  chosen 
to  serve  as  a  rigorous  test  of  the  calculation. 
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Fig.  5.32.  Comparison  of  Albedo  Calculations  by 
Simon-Clifford  Technique  with  Thermal-Neutron  Fluxes 
Meosured  in  Straight  Cylindrical  Steel-Walled  Ducts 
Through  a  Water  Shield:  Plane  Source  of  Thermal  Neu¬ 
trons.  (From  Horton  and  Halliday,  ref.  39.) 


ORNL-DWG  66-12477 


TOP  BUNKER 


Fig.  5.33a.  Schematic  of  Two  Concrete  Bunkers  (12- 
ft  Cubes)  Connected  by  3- X  8-ft  Tunnel. 
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since  the  fluxes  and  doses  for  detector  locations 
greater  than  a  few  feet  from  the  duct  mouth  would 
then  be  due  entirely  to  wall-scattered  radiation. 
The  measured  absolute  energy  spectrum  of  the 
beam  was  used  for  the  source  spectrum  in  the 
calculations. 

Figure  5.34  shows  a  comparison  of  the  calcu¬ 
lated  and  measured  results  for  fast  neutrons.  In 


general,  the  agreement  averages  better  than  20% 
through  five  orders  of  fast-neutron  dose  atten¬ 
uation.  In  the  case  of  subcadmium  or  thermal 
fluxes  that  result  from  an  incident  beam  of  sub¬ 
cadmium  neutrons  (Fig.  5.35),  the  agreement  is 
within  20%  in  the  second  and  third  legs  and  within 
5%  for  most  of  the  straight  ducts. 
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Fig.  5.336.  Comparison  of  Monte  Carlo  Albedo  Calculations  with  Measurements  of  Thermal-Neutron  Fluxes 
in  Three-Legged  Rectangular  Concrete  Tunnel.  (From  Coin,  ref.  17.) 
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In  order  to  test  the  accuracy  of  the  slowing- 
down  model  used  in  the  code,  calculations  were 
also  made  of  the  thermal-neutron  flux  resulting 
from  all  source  neutrons  with  energies  above  cad¬ 
mium  cutoff.  The  comparison  between  calculations 
and  experiment  shown  in  Fig.  5.36  indicates  that 
the  agreement  is  about  the  same  as  that  obtained 
for  incident  subcadmium  neutrons.  The  statistical 
errors,  number  of  wall  backscatterings,  and  run¬ 
ning  times  for  the  same  number  of  source  his¬ 
tories  were  comparable  in  the  two  calculations. 
It  should  be  noted  that  the  calculations  were  nor¬ 
malized  in  each  case  to  the  measured  number  of 
neutrons  incident  on  the  wall  of  the  duct  from 
the  source  beam. 

Calculations  were  also  made  of  the  epicadmium 
or  nonthermal  spectra  at  several  locations  along 
the  center  line  of  the  three-legged  duct,  and  these 

Fig.  5,34.  Comparison  of  AMC  Albedo  Calculations 
with  Measurements  of  Fast-Neutron  Dose  Rates  in  One- 
and  Two-Legged  3-ft-square  Concrete  Ducts.  The 
source  in  the  experiment  was  a  collimated  beam  from 
the  Tower  Shielding  Reactor  II,  (From  Maerker,  ref. 

14.) 
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Fig.  5.35.  Comparison  of  AMC  Albedo  Calculations  with  Measurements  of  Thermal-Neutron  Fluxes  in  One-, 
Two-,  and  Three-Legged  3-ft-square  Concrete  Ducts  Doe  to  Incident  Subcadmium  Neutrons.  (From  Maerker,  ref.  14.) 
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1.45  X  10“®(1.2  X  10“®) 
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The  dose  rates  in  parentheses  are  those  obtained  when  the  center  of  detection  was  assumed  to  be  at  the  center 
of  the  polyethylene  ball  rather  than  at  its  leading  edge. 
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spectra,  together  with  the  flux-to-multicollision 
dose  factors  of  Snyder  and  Neufeld,^^  were  used 
to  calculate  the  epicadmium  multicollision  dose 
rates  at  the  same  locations.  The  results  are 
compared  in  Table  5.5  with  measurements  made 
with  a  dosimeter  designed  to  give  a  response 
proportional  to  the  Snyder-Neufeld  multicollision 
dose  curve  over  the  range  of  interest.  The  do¬ 
simeter  was  a  spherical  BF^  counter  placed  in¬ 
side  a  10-  or  12-in. -diam  polyethylene  ball  covered 
with  cadmium. 

5.1.5.  EMPIRICAL  CORRELATION 
FOR  RECTANGULAR  DUCTS 

From  the  background  of  experimental  and  an¬ 
alytical  data  available,  Ingold  and  Huddleston^* 
developed  an  empirical  formula  for  estimating 
the  attenuation  of  gamma  rays  in  rectangular  ducts. 
By  best  fit  to  the  data,  the  formula  was  found  to  be 

D  W'2.864 

—  =0.214 - ,  (5.51) 

D„  T  2.534  j  2.667  p0.710 

“  ^1  ^2  ^0 

where 

D  =  dose  rate  on  axis  of  duct  at  point  of  interest, 
Dg  =  dose  rate  1  ft  from  the  source, 


H  =  height  of  duct  (ft), 

H'  =  width  of  duct  (ft), 

Lj  =  length  of  first  leg  measured  from  the  source 
to  the  center  of  the  first  corner  (ft), 

=  length  of  second  leg  measured  from  the 
center  of  the  corner  to  the  detector  (ft). 

Eg  =  gamma-ray  source  energy  (MeV). 

Use  of  this  formula  is  limited  to  problems  which 
meet  the  following  criteria: 

Point  source  at  duct  entrance,  0.662  MeV  = 
Eg  S  6.0  MeV 

1.0^  6.0  ft 

1.0^  6.0  ft 

2  ^  ^  36  ft 

1  ^  H/W  ^  2 
L/E^6 
L^/hS6 
2 

L^/W^2 

While  this  formula  is  useful  for  the  case  of  a 
point  source,  it  should  be  applied  with  caution 
to  the  case  of  a  broad-beam  source. 

Results  obtained  with  Eq.  5.51  are  shown  in 
Figs.  5.21  and  5.22  (Section  5.1.4.). 


5.2.  Methods  for  Calculating  Radiation  Transmission 
Through  Shields  Containing  Voids 


5.2.1.  RAY  ANALYSIS  TECHNIQUE  Such  irregularities,  which  can  vary  in  size  and 

FOR  SINGLE  VOIDS  shape,  are  usually  referred  to  as  “voids.”  In 

general,  voids  do  not  extend  either  to  the  source 
There  will  be  situations  in  which  a  shield  will  or  to  the  shield  surface,  but  methods  applied  to 

contain  irregularities  that  are  not  classifiable  as  them  can  also  be  used  to  treat  depressions  in 

ducts  but  which,  like  ducts,  cause  the  flux  at  shield  surfaces. 

nearby  points  on  the  shield  surface  to  be  increased  The  technique  most  frequently  applied  to  single 
simply  because  the  amount  of  attenuation  offered  voids  is  the  ray  analysis  technique,  which  is 

by  that  region  of  the  shield  has  been  reduced.  described  in  Section  5.1.2  in  connection  with  duct 
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transmission  problems.  As  is  pointed  out  there, 
the  basic  assumption  of  the  ray  analysis  technique 
is  that  the  radiation  transmission  is  a  function  only 
of  the  path  lengths  through  each  material  or  void 
encountered  along  a  straight  line  between  the 
source  point  and  the  point  of  interest.  When  ap¬ 
plied  to  voids,  the  increase  in  the  flux  on  the 
shield  surface  is  determined  by  the  difference 
in  the  fluxes  calculated  with  and  without  the  void 
present. 

Consider,  for  example,  the  case  of  a  disk-shaped 
void  in  a  semi-infinite  slab  shield  adjacent  to  an 
infinite  plane  gamma-ray  source  with  the  ends  of 
the  disk  parallel  to  the  infinite  dimension  as  is 
shown  in  Fig.  5.37.  Let  ^(L)  be  the  flux  of  energy 
jE  at  a  point  on  the  surface  of  the  shield  with  no 
void,  and  0(L  —  t)  be  the  flux  with  the  void  present. 
For  the  latter  the  point  of  interest  is  the  intersec¬ 
tion  of  the  center  line  of  the  void  with  the  shield 
surface.  The  influence  of  the  void  is  estimated  by 
performing  a  point  kernel  integration  over  the 
source  plane  to  obtain  <I>(L)  and  $(£  —  t). 

To  calculate  the  flux  for  the  plane  shield  of 

thickness  L,  define  =  2  over  L  normally 

i 

through  the  shield.  The  contribution  of  a  ring- 
shaped  source  with  a  particle  emission  rate  of  N ^ 
at  an  angle  6  from  the  detector  would  then  be 


^ 0  dr  —b  seed 

= - e  ‘ 

2  r 


(5.52) 
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Fig.  5.37.  Geometry  for  Calculating  Effect  of  Disk- 
Shaped  Void  on  Radiation  Transmission  Through  Shield. 


Integrating  over  the  plane  source  out  to  an  angle 
6  gives 


A^o 

$(0)=_  /  - dt  ,  (5.53) 

2  t 

1 


where  t  =  sec  6.  Since  by  definition 


(5.54) 


then 


^(6)  =-^  [£j(6j)  -  E,(bi  sec  0)]  ,  (5.55) 

and  integrating  over  the  entire  source  plane  gives 

<D(L)  =_  Ej(b^)  .  (5.56) 

If  the  gamma-ray  number  buildup  through  is  B 
then  0(L)  may  be  estimated  by 

^(L)  =^-  E^ib^)  .  (5.57) 

The  buildup  factor  used  here  as  a  simple  multi¬ 
plier  is  only  symbolic  since  the  buildup  factor  is 
really  part  of  the  kernel  that  is  integrated  over 
all  space.  As  is  pointed  out  in  Chapter  3,  when 
the  Taylor  form  of  the  buildup  factor  is  used,  the 
form  of  these  equations  is  unchanged  and  the 
buildup  effect  is  simply  included  by  modifying 
the  arguments  of  the  exponential  functions.  Actu¬ 
ally  in  these  problems  it  is  sufficiently  accurate 
to  use  the  equation  as  it  stands  and  to  determine 
only  on  the  basis  of  the  slant  penetration  (Q  to 
P)  that  grazes  the  void.  For  B an  average  of 
B  j  and  of  the  buildup  factor  based  on  the  minimum 
path  length  normal  to  the  void  is  usually  quite 
adequate.  Using  the  arithmetic  mean  will  always 
produce  a  conservative  result  since  this  gives 
equal  weight  to  the  buildup  factor  from  a  region  of 
lesser  importance. 

Now  to  calculate  the  flux  for  a  shield  in  which 
a  disk-shaped  void  intercepts  at  an  angle  d  from 
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the  point  P,  define  =  2  over  (L  —  t).  Then 
for  angles  less  than  d 


<D  = - 

^  2 


[E^(b^)-E^(b^sece)]  ,  (5.58) 


and  for  angles  greater  than  0 

= — : — 


2  2 
Thus  the  total  flux  is 


(5.59) 


$(L  -  f)  =  'l>i  +  'I'a 


£j(hj  sec  d) 


+  [E  j(hp  -  E  ^(b^  sec  d)]  ,  (5.60) 

and  the  difference,  ^(L  —  t)  —  $(L),  is  the  in¬ 
crease  in  the  flux  at  point  P  due  to  the  presence 
of  the  void.  Because  P  is  centered  over  the  void, 
it  will  be  the  point  of  maximum  flux  increase. 

For  points  offset  from  the  center,  for  example, 
for  point  P'  in  Fig.  5.37,  Eq.  5.60  will  not  be 
valid.  If  the  source  is  isotropic,  the  increase  in 
the  contribution  at  P  '  by  a  source  surface  increment 
dS  caused  by  the  void  is  given  by 

N^dS 

^i  = - : 

27t(L  sec  j8)^ 

X  [S3  g-^LsecS]  ^  (5  61) 

where  B  ^  and  B  ^  are  the  buildup  factors  along  the 
path  dS  ^  P'  with  and  without  the  void  present 
respectively.  If  the  source  is  a  cosine  source, 
the  increase  is  given  by 

N  dS  cos  ^ 

^c  =  - - : 

77‘(L  sec  fSy 

X  _  (5  62) 

The  integral  of  Eq.  5.61  or  5.62  over  the  appropri¬ 
ate  source  area  will  give  the  total  increase  in  flux 
at  P'  due  to  the  presence  of  the  void.  The  appro¬ 
priate  source  area  includes  all  points  for  which  a 


ray  to  P '  passes  through  the  void.  In  most  prac¬ 
tical  situations  this  integral  must  be  evaluated 
numerically.  The  value  of  the  integral  at  repre¬ 
sentative  points  will  serve  to  determine  the  shape 
and  thickness  of  the  shield  patch  required  to 
negate  the  effect  of  the  void. 

A  simple  approximate  form  of  ray  analysis  given 
by  Tonks'*^  may  be  applied  to  certain  problems 
involving  voids  or  depressions  in  a  very  thick 
shield.  This  method  considers  first  the  flux  at  a 
point  on  the  surface  of  a  uniform  infinite  slab 
shield  whose  thickness  is  equal  to  the  ray  of 
minimum  length  through  a  void  or  depression  in 
the  shield.  It  then  accounts  for  the  reduction  in 
flux  due  to  the  longer  material  path  of  nearby  rays 
terminating  at  the  same  point.  The  method  is 
based  on  the  fact  that  the  formula  for  the  flux 
through  an  infinite  slab, 

<h=_E^(^0-  (5.63) 

may  be  approximated  by 
N 

(D  =  — (5.64) 

2fit 

for  fit  »  1  with  an  error  that  is  approximately  equal 
to  l/flt. 

The  application  of  Tonks’  approximation  is  as 
follows.  Consider  a  slab  shield  with  a  void  such 
as  that  shown  in  Fig.  5.38.  Let  the  ray  passing 
through  Q  and  P  be  the  path  of  minimum  material 
penetration,  and  let  it  coincide  with  the  z 

axis.  In  many  cases  nearby  rays  may  be  expressed 
as 


f  ^  ^  ^  ’  (5-65) 

where  6  and  6  are  the  angles  between  the  mini- 
mum  ray  and  a  nearby  ray  projected  onto  the  xz 
and  yz  planes  respectively  and  a  and  b  are  deter¬ 
mined  by  the  specific  geometry.  Substituting  Eq. 
5.65  into  an  approximation  analogous  to  Eq.  5.64 
yields  for  the  flux  <1) '  at  P 


=/V„ 


II  ’ 


(5.66) 
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Other  expansions  of  i 

derived  which  would  lead  to  equations  of  the 
approximate  form  of  Eq.  5.66. 


Fig.  5.38.  Geometry  for  Calculating  the  Effect  of  a 
Void  in  a  Shield  by  Tonhs'  Approximation. 


if  both  a  and  b  are  greater  than  zero.  The  ratio  of 
the  flux  estimated  by  Eq.  5.66  to  that  estimated 
by  Eq.  5.64  (for  a  slab  of  minimum  penetration) 
is 


<&  \/a6  ’ 


(5.67) 


5.2.2.  FLUX  PERTURBATION  TECHNIQUE 
FOR  SINGLE  VOIDS 

Routs'*^  treated  the  problem  of  voids  on  the 
basis  of  one-velocity  transport  theory.  Instead  of 
using  the  boundary  conditions  of  the  differential 
transport  equation  to  describe  the  void,  he  ac¬ 
counted  for  it  by  the  way  in  which  the  kernel  of 
the  integral  equation  is  defined.  The  effect  of  the 
void  is  then  calculated  as  a  perturbation  of  the 
flux  in  a  homogeneous  medium  (no-void  case).  The 
method  is  applicable  to  either  gamma  rays  or 
neutrons. 

Using  this  approximation  to  the  transport  equa¬ 
tion,  Routs  investigated  the  effect  of  spherical 
voids  in  water  and  obtained  good  agreement  with 
experiment  when  he  chose  certain  constants  in 
the  solution  by  fitting  to  the  experimental  data  at 
one  point.  The  expression  derived  for  a  spherical 
void  in  the  geometry  shown  in  Fig.  5.39  has  the 
form 


which,  provided  that  >  1,  gives  the  increase 

in  effectiveness  at  P  of  the  shield  with  a  void 
over  that  of  a  simple  slab  whose  thickness  is 
equivalent  to  the  minimum  ray. 

The  limit  of  6^  and  0^  for  which  Tonks’  approx¬ 
imation  is  valid  is  set  hy  0  —  6  <  cl,  where  a 

_  A  y 

is  given  by  a  =  y/2/fit^.  As  an  example,  for  a 
spherical  void  of  radius  whose  center  lies  at 
a  depth  R  from  the  outside  surface  of  a  slab  shield 
of  thickness  T,  it  can  be  shown  that  the  length  of 
a  ray  near  the  minimum  ray  is  given  by 

/r  R^\  ,  ^0  , 

t^t^  +  {-+  —  )0l  +  -0l,  (5.68) 

“  \2  R^.J  2  y 

and  the  increased  effectiveness  at  the  minimum-ray 
terminus  is  given  by 


O' 


t 


0 


0 


'T  R^\  ^0 


1/2 


(5.69) 


0(P) 


where 


(5.70) 


0(P)  =  flux  at  point  P, 

B  =  constant  in  the  approximation  to  the 
angular  distribution  of  the  scattered 
flux  used  by  Routs,  that  is. 


/  dfl  0(n)2p(n'  -fi) 


=  B  cos  0, 


o<d: 


77 

=  0,  —  <  0  <  77, 

2 

2  =  total  macroscopic  cross  section, 

p(Q'  “*0)  =  probability  of  scattering  from  direc¬ 
tion  fl  '  to  direction  0, , 

V  =  parameter  in  the  expression  for  flux 
distribution  through  the  shield  in  the 
absence  of  voids  given  by  = 

F<I'(0)  e-3/*. 
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R  =  radius  of  void, 

L  =  shield  thickness, 
b  =  number  of  relaxation  lengths  in  the 
shield  material, 

=  cos  6. 


A  best  value  of  k  =  4.5  was  chosen  by  Kouts  by 
fitting  experimental  data  for  gamma  rays  at  2R/L  = 

1. 

Results  obtained  with  this  equation  are  compared 
in  Fig.  5.40  with  experimental  data  for  gamma  rays. 
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Fig.  5.39.  Geometry  for  Calculating  Effect  of  Void  in 
Shield  by  Kouts'  Technique. 
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Fig.  5.40.  Comparison  of  Flux  Increase  Due  to  Spheri¬ 
cal  Void  in  Water  Shield  Calculated  by  Kouts'  Technique 
with  That  Determined  by  Experiment.  (From  Kouts, 
ref.  43.) 


5.2.3.  TECHNIQUE  FOR  SMALL  RANDOMLY 
SPACED  VOIDS 


For  the  case  of  a  large  shield  containing  small 
randomly  spaced  voids  which  are  small  in  distance 
across  relative  to  the  relaxation  length  of  the 
radiation  (for  example,  a  gamma-ray  shield  con¬ 
sisting  entirely  of  iron  shot),  the  effect  of  the 
voids  on  the  radiation  transmission  is  obtained 
simply  by  increasing  the  relaxation  length  used 
in  the  calculation  of  the  shield  attenuation  by  the 
ratio  of  the  average  density  of  the  shield  to  the 
density  of  the  solid  material  used  in  the  shield. 

When  the  voids  are  not  negligibly  small,  however, 
it  can  be  shown  that  on  a  statistical  basis  the 
penetrating  radiation  will  be  greater  than  would 
be  predicted  by  using  an  average  density.  Con¬ 
sider,  for  example,  two  normal  rays  through  the 
shield,  one  which  is  the  average  penetration  ( 
plus  an  increment  §  and  one  which  is  shorter  than 
average  by  the  same  increment.  Their  average 
penetration  is 

1  =  e-^  '  cosh  Su  >  . 

2 

Thus,  although  their  average  penetration  distance 
is  t,  their  average  penetration  is  greater  than 

This  problem  has  been  treated  statistically  by 
Smith, who  arrived  at  an  effective  thickness  t 

e 

given  by 


f(l  -  V) 

-  f 

1  +  0.55sv^ 


(5.71) 


where  t  is  the  actual  thickness,  v  is  the  fraction 
of  the  volume  which  is  void,  and  s  is  the  average 
distance  between  voids  measured  in  relaxation 
lengths.  An  equivalent  expression  for  an  effective 


attenuation  coefficient  is 
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or,  generalized  to  the  case  where  the  voids  contain 
a  material  of  attenuation  coefficient  ,  is 


/r(l  -  v) 

1  +  0.55  (s/^) 


(5.72) 


ii'  +  a  -  v)iti  -  lii 
1  +  0.55  (s/,x)  ill -111 


(5.73) 


Appendix  5A.  Machine  Programs 


Several  machine  programs  that  utilize  either 
the  ray  analysis  or  the  analog  Monte  Carlo  tech¬ 
nique  described  in  Sections  5.1.2  and  5.1.3  re¬ 
spectively  are  broadly  outlined  below  to  aid  the 
reader  in  deciding  which  of  these  particular  pro¬ 
grams  most  nearly  applies  to  his  particular  sit¬ 
uation.  It  is  emphasized  that  other  good  computer 
codes  are  available  and  new  ones  are  constantly 
being  developed.  Further  information  on  shielding 
programs,  including  the  availability  of  the  pro¬ 
gram  decks,  libraries,  and  utilization  instructions, 
may  be  obtained  from  the  USAEC  Radiation 
Shielding  Information  Center  at  Oak  Ridge  National 
Laboratory,  P.  O.  Box  X,  Oak  Ridge,  Tennessee. 


5A.1.  RAY  ANALYSIS  PROGRAMS 

Some  common  features  of  computer  programs 
utilizing  the  ray  analysis  or  point-kernel  tech¬ 
niques  are  as  follows: 

(1)  The  source-shield  geometry  is  described  as 
a  series  of  simple  geometric  regions. 

(2)  The  source  is  divided  into  a  number  of  regions, 
and  the  radiation  emitted  from  each  region 
is  assumed  to  come  from  a  point  concentrated 
at  its  center. 

(3)  The  flux  or  some  response  to  flux  (exposure 
dose,  heating,  etc.)  is  calculated  at  designated 
points  by  summing  the  contribution  from  each 
of  the  source  points. 


C-17  and  L-63  Codes.  —  The  General  Dynamics 
C-17  and  L-63  shield  penetration  programs'*®  are 
identical  except  for  the  geometry  subprogram.  Ac¬ 
ceptable  geometries  for  the  C-17  program  are 
limited  to  frusta  of  rectangular  pyramids  and 
coaxial  cylinders  and  their  annuli.  The  L-63  pro¬ 
gram  is  more  inclusive  in  that  it  accepts  cylinders, 
and  their  annuli,  that  are  defined  about  arbitrary 
axes,  sectors  of  these  cylinders,  and  frusta  of 
pyramids  whose  bases  are  quadrilaterals,  spheres, 
hollow  spheres,  hemispheres  and  spherical  sec¬ 
tors.  The  L-63  program  also  accepts  regions 
within  regions  and  regions  within  regions  within 
regions  in  which  the  geometry  types  can  be  varied. 

The  geometry  routine  determines  the  intercepted 
distance  in  each  material  along  each  source- 
detector  ray  by  an  iterative  stepping  point.  These 
distances  are  used  in  the  gamma-ray  and/or  neu¬ 
tron  routines  for  calculation  of  the  attenuation 
along  each  ray. 

The  gamma-ray  routine  calculates  the  differential 
gamma-ray  energy  spectrum  at  a  point  detector 
for  each  point  source.  The  uncollided  flux  cal¬ 
culation  is  carried  out  in  a  straightforward  manner. 
The  scattered  flux  is  evaluated  by  means  of  the 
NDA  moments  method  data,'*®  which  give  the 
spectra  of  scattered  gamma  rays  due  to  mono- 
energetic  point  isotropic  sources  in  an  infinite 
medium  as  a  function  of  penetration  depth.  In 
applying  these  data  to  finite  media,  the  influence 
of  material  interfaces  must  be  discounted.  The 
total  spectrum  calculated  at  a  detector  point  is 
the  sum  of  the  spectra  from  the  individual  source 
points.  The  total  spectrum  may  then  be  multiplied 
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point-wise  by  either  flux-to-exposure-dose  or  flux- 
to-heat  conversion  factors.  For  materials  which 
are  mixtures  of  elements,  an  effective  atomic 
weight  is  determined  and  the  differential  number 
spectra  are  found  by  curve  fitting  to  the  data  of 
ref.  46.  A  formula  for  calculating  effective  atomic 
numbers  is  given  in  the  utilization  manual.'*® 

Basic  data  on  neutron  penetration  are  not  nearly 
so  comprehensive  as  the  data  on  gamma-ray  at¬ 
tenuation.  In  order  to  determine  the  neutron  spec¬ 
trum  from  an  individual  point  source,  a  reference 
material  must  be  selected  for  which  basic  pene¬ 
tration  calculations  have  been  performed.  All 
materials  along  the  ray  are  then  assumed  to  be 
the  reference  material.  The  penetration  distances 
for  each  material  are  weighted  in  accordance 
with  the  ratio  of  the  material  removal  cross  section. 

QAD  Program.  —  The  Los  Alamos  Scientific 
Laboratory  program  QAD'*’^  accepts  any  geometry 
which  can  be  described  by  a  set  of  cartesian 
quadratic  equations.  It  cannot  calculate  differ¬ 
ential  gamma-ray  spectra.  A  single  material 
buildup  factor  is  used  in  calculating  the  dose 
at  a  point  due  to  scattered  gamma  rays.  Neutron 
attenuation  is  determined  by  interpolation  of  basic 
penetration  information  for  a  single  base  material 
in  a  manner  similar  to  the  General  Dynamics  C-17 
program  (see  above). 

GE  04-4  Program.  —  The  General  Electric  Pro¬ 
gram  04-4^®  may  be  used  to  calculate  neutron 
and  gamma-ray  dose  rates  around  complex  re¬ 
actor-shield  assemblies.  It  will  not  compute 
energy  spectra  for  either  type  of  radiation.  Ac¬ 
ceptable  geometries  are  those  which  may  be  de¬ 
scribed  by  rotating  rectangles  or  trapezoids  about 
an  axis  or  by  translating  rectangles. 

Scattered  gamma-ray  dose  is  determined  by  a 
fit  to  buildup  factors  taken  from  the  NDA  moments 
method  data."*®  The  buildup  as  a  function  of 
penetration  is  approximated*®  by  a  series  of  poly¬ 
nomial  coefficients  for  each  material  and  energy 
group.  Buildup  factors  for  laminated  shields  are 
computed  using  the  formulas  derived  for  lead-water 
and  iron-water  combinations  by  Kalos.®®  The 
assumption  is  made  that  these  expressions  may 
be  applied  to  other  material  combinations. 

Fast-neutron  attenuation  is  calculated  using 
an  approach  suggested  by  Albert  and  Welton.  ®  * 
In  this  theory  collisions  with  hydrogen  are  treated 
as  absorptions,  and  an  approximate  relation  is 
used  to  describe  the  total  cross  section  as  a 


function  of  energy.  For  the  calculation  of  at¬ 
tenuation  by  materials  other  than  hydrogen,  ex¬ 
ponential  attenuation  is  assumed  and  the  cross 
sections  are  treated  as  energy-independent  ad¬ 
justable  parameters  to  be  determined  by  best  fit 
to  experimental  data.  This  treatment  is  based 
on  the  assumption  that  heavy  materials  along 
the  ray  are  followed  by  sufficient  hydrogen  to 
make  removal  theory  valid. 

5A.2.  ANALOG  MONTE  CARLO  PROGRAMS 

LOS  Code.  —  The  LOS  Monte  Carlo  computer 
routine®^  may  be  used  for  analysis  of  radiation 
streaming  down  cylindrical  ducts.  The  atomic 
density  and  cross  sections  of  the  elements  in 
the  material  around  the  duct  are  used  as  input, 
and  the  duct  is  defined  as  a  region  of  zero  or  re¬ 
duced  density,  as  appropriate. 

The  restriction  on  geometry  is  that  all  surfaces 
must  be  cylindrical  and  coaxial  between  any  two 
planes.  Since  geometry-dividing  planes  need  not 
be  parallel  and  the  axis  may  vary  freely  between 
consecutive  planes,  multibend  ducts  of  consid¬ 
erable  complexity  may  be  defined.  Material  com¬ 
position  may  vary  radially  or  axially  with  regions 
of  constant  composition  being  separated  by  cy¬ 
lindrical  surfaces  or  planes. 

Some  of  the  more  sophisticated  features  of  the 
LOS  procedure  include:  biased  sampling  from 
source  angular  distributions,  variance  reduction 
by  splitting  and  Russian  roulette  at  boundaries, 
and  the  ability  to  calculate  flux  or  dose  rates 
at  nondirectional  point  detectors  by  means  of 
statistical  estimation. 

Allowable  neutron  interactions  include  any  com¬ 
bination  of  isotropic  or  anisotropic  elastic  and 
inelastic  scattering  in  the  laboratory  or  the  center- 
of-mass  system.  All  significant  photon  inter¬ 
actions  are  allowed. 

ADONIS  Code.  -  The  ADONIS  Monte  Carlo  pro¬ 
gram  *  *  may  be  used  to  analyze  neutron  or  gamma- 
ray  transport  in  regions  which  may  be  expressed 
as  rectangular  parallelepipeds.  Either  straight-  or 
bent-duct  streaming  may  be  calculated  with  this 
code.  The  total  shield  and  duct  geometry  may 
be  described  by  using  up  to  80  rectangular  paral- 
lelepipedal  regions.  Only  right-angle  bends  are 
allowed,  although  bends  at  other  angles  may  be 
simulated  by  a  series  of  short  rectangular  regions 
provided  that  the  80-region  limit  is  not  exceeded. 
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Like  the  LOS  code,  this  program  also  employs 
splitting  and  Russian  roulette.  Flux  is  calculated 
on  the  basis  of  track  length  per  unit  volume  in 
specified  regions.  Flux  may  be  multiplied  by 
any  response  function  in  a  region  to  obtain  dose, 
heating  rate,  interaction  rate,  or  secondary  source 
strength. 

A  later  version  of  the  ADONIS  code,  called 
UNC-SAM,  has  been  used  by  United  Nuclear  to 
determine  neutron  streaming  in  a  passageway  due 
to  a  reactor  source.  Results  obtained  with  this 
version  are  reported  to  be  in  excellent  agreement 
with  experiments. 

COHORT  Code.  -  The  COHORT  Monte  Carlo 
code®'*  was  written  for  the  analysis  of  neutrons 
and  gamma  rays  in  complex  geometries.  Primary 
advantages  of  the  method  are  the  ability  to  cal¬ 
culate  flux  in  volume  or  point  detectors  and  the 
fact  that  secondary  sources  may  be  analyzed 
without  resort  to  external  hand  calculation.  The 
program  accepts  any  geometry  which  may  be  de¬ 
scribed  by  bounding  spherical,  hyperbolic,  conical, 
cylindrical,  or  plane  surfaces.  All  surfaces  except 
planes  must  be  symmetric  about  a  line  parallel 
to  the  Z  axis.  COHORT  is  composed  of  a  family 
of  six  codes;  SOI  for  generating  source  particles, 
S02  for  generating  secondary  source  particles, 
HOI  for  generating  collision  histories,  AOl  for 
analyzing  histories  to  determine  flux  at  a  point, 
A02  for  analyzing  histories  to  determine  track 
length  per  unit  volume,  and  COl,  an  auxiliary 
program,  which  is  used  to  interpret  and  print  out 
binary  tapes  generated  by  other  programs.  All 
codes  are  written  in  FORTRAN  IV.  Exponential 


transformation  and  splitting  options  are  available 
for  variance  reduction  in  the  HOI  routine. 

05R  Code.  -  The  Oak  Ridge  National  Labora¬ 
tory  05R  Monte  Carlo  code®®  is  a  versatile  pro¬ 
gram  for  the  analysis  of  neutron  transport.  It 
provides  a  particularly  detailed  treatment  of  cross 
sections.  Geometry  description  is  limited  only 
by  the  requirement  that  surface  boundaries  be 
represented  by  a  general  quadratic  or  a  plane. 
Up  to  16  media  may  be  represented.  Component 
parts  of  the  program  include  the  following: 

XSECT,  which  consists  of  a  package  of  nine 
basic  cross-section  handling  codes.  Various 
components  of  this  package  and  a  separate 
program  called  LEGENDRE  are  used  in  the 
preparation  and  modification  of  a  master  cross- 
section  tape  and  two  tapes  used  by  05R  to  ob¬ 
tain  cross  section  and  angular  distribution  data. 

GEOM,  which  is  the  geometry  subroutine. 

05R,  the  “heart”  of  the  program,  which  is  used 
to  generate  collision  parameters  which  are  stored 
on  a  data  tape. 

The  program  is  incomplete  in  that  the  user  must 
provide  his  own  source  generation  and  analysis 
routines.  Additionally,  the  user  must  provide  a 
program  for  processing  inelastic  scattering  events. 
The  user  may  provide  any  desired  variance  re¬ 
duction  in  the  analysis  routine.  Splitting  and 
Russian  roulette  are  allowed  at  boundaries  which 
may  be  arbitrarily  assigned  in  the  geometry  routine. 
Unit  sources  (all  starting  parameters  the  same) 
may  be  analyzed  without  use  of  a  separate  source 
routine. 
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